Integrable Systems and Factorization Problems 

M. A. Semenov-Tian-Shansky 



The present lectures were prepared for the Faro International Summer School on 
Factorization and Integrable Systems in September 2000. They were intended for par- 
ticipants with the background in Analysis and Operator Theory but without special 
knowledge of Geometry and Lie Groups. The text below represents a sort of compro- 
mise: it is certainly impossible not to speak about Lie algebras and Lie groups at all; 
however, in order to make the main ideas reasonably clear, I tried to use only matrix 
algebras such as and its natural subalgebras; Lie groups used are either GL(n) and 
its subgroups, or loop groups consisting of matrix-valued functions on the circle (pos- 
sibly admitting an extension to parts of the Riemann sphere). I hope this makes the 
environment sufficiently easy to live in for an analyst. The main goal is to explain how 
the factorization problems (typically, the matrix Riemann problem) generate the entire 
small world of Integrable Systems along with the geometry of the phase space, Hamilto- 
nian structure, Lax representations, integrals of motion and explicit solutions. The key 
tool will be the classical r-matrix (an object whose other guise is the well-known Hilbert 
transform). I do not give technical details, unless they may be exposed in a few lines; on 
the other hand, all motivations are given in full scale whenever possible. I hope that this 
choice agrees with the spirit of the Faro School and will help to bridge the gap between 
different branches of Mathematical Analysis. 



1. Introduction 



The story of the discovery of the modern theory of Integrable Systems is cer- 
tainly too long (and too well-known), and I can hardly add anything new; so let 
me tell just a few words before addressing the bulk of the subject. The study 
of completely integrable systems goes back to the classical papers of Euler, La- 
grange, Jacobi, Liouville and others on analytical mechanics. By the end of the 
XlX-th century all interesting examples seemed to have been exhausted, and the 
interest has shifted to the qualitative study of chaotic behaviour. The new age 
in the study of integrable systems has begun with the famous paper [1] on the 
KdV equation, which was the first example of an infinite dimensional dynamical 
system with nontrivial but highly regular behaviour and with a rich excitation 
spectrum. Amplifying the earlier results of Gardner, Greene, Kruskal and Miura 
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and of Peter Lax [2] Faddeev and Zakharov [3] have shown in 1971, exactly 30 
years ago, that the KdV equation is in fact a completely integrable Hamiltonian 
system in a technical sense. Within a short time, many more examples have 
been discovered, notably, the sine-Gordon equation, the first ever example of a 
relativistic completely integrable system [4]. These discoveries were particularly 
exciting in view of the possible physical applications: while it was of course clear 
that "generic" nonlinear equations are non-intcgrable, it has been argued that 
Fundamental Physics always deals with highly non-generic c;quations which might 
be integrable in some sense or other. These initial hopes have been fulfilled only 
partially; one major obstacle is that the new technique does not apply to (nat- 
ural) non-linear equations in realistic space-time dimension. On the other hand, 
the mathematics of complete integrability has proved to be extraordinarily rich, 
bringing together Functional Analysis, Algebraic Geometry, Lie Groups, Repre- 
sentation Theory, Symplectic Geometry and much more. The analytic machinery 
used in the initial papers was based on the Inverse Scattering Problem; the subse- 
quent developments allowed to single out the basic geometric ideas of the theory 
and to provide a unified basis for different examples. One of the key ingredients of 
this geometric approach are Infinite Dimensional Lie Groups; in some loose sense, 
one can say that Integrable Systems always possess some rich hidden symmetry. 
One may recall that one of the original motivations of S.Lie has been the use of 
continuous transformation groups for the study of differential equations. With the 
modern methods at hand, we came closer to that goal; it has now become clear, 
however, that almost all classical examples of integrable mechanical systems from 
the XIX century textbooks, as well as the infinite dimensional systems associated 
with integrable PDE's are related not to finite-dimensional Lie groups, but rather 
to their infinite- dimensional analogs. 

The general geometric construction that we shall discuss below allows to unify 
the following characteristic features that are typical for all examples known so far: 

(i) The equations of motion are compatibility conditions for a certain auxiliary 
system of linear equations. 

(ii) They are Hamiltonian with respect to a natural Poisson bracket. 

(iii) Integrals of motion are spectral invariants of the auxiliary linear operator. 
They are in involution with respect to the Poisson bracket referred to above. 

(iv) The solution of the equations of motion reduces to some version of the 
Riemann-Hilbert problem. 

Depending on the nature of the auxiliary linear problem, the associated nonlinear 
equations may be divided into the following three groups: 

a) Finite-dimensional systems, 

b) Infinite-dimensional systems with one or two spatial variables, 

c) Integrable systems on one-dimensional lattices. 
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In case (a) the auxiliary linear problem is the eigenvalue problem for a finite- 
dimensional matrix (possibly depending on an additional parameter). In case (b) 
the associated linear operator is differential. In case (c) it is a difference operator. 

As it happens, the key properties 1-4 referred to above are corollaries of a single 
general theorem which may be adapted to numerous concrete applications. The 
original idea of this theorem is due to B.Kostant and M.Adler; its important am- 
plification brings in the notion of classical r-matrix which provides a link between 
abstract Riemann-Hilbert problems and the ideas of Quantum Group Theory. The 
statement and proof of this theorem are particularly simple for systems of types (a) 
and (b). (Lattice systems require a special treatment, since the associated Poisson 
brackets belong to a different and more sophisticated class. We shall discuss this 
case later on, but it is natural to begin with the simpler cases (a) and (b)). 

One more word of caution: while it is aesthetically very attractive to deduce 
a large variety of examples together with their explicit solutions from a simple 
general construction, there is one important disadvantage: for a given dynamical 
system (even if it is known to be completely integrable!) it is very difficult to 
tell a priori, what is the underlying Lie group, or Lie algebra. The practical way 
around this difficulty is to look at various examples associated with different Lie 
algebras; with some skill, one manages to recognize among these examples both 
classical and new integrable systems which admit physical interpretation. The list 
of interesting Lie algebras includes: 

1. Finite dimensional semisimple Lie algebras. Associated integrable systems 
include open Toda lattice and other finite dimensional Hamiltonian systems 
which may be integrated in elementary functions (rational functions of expt, 
or t, where t is the time variable). 

2. Loop algebras, or affine Lie algebras. Associated integrable systems are 
finite dimensional Hamiltonian systems which may be integrated in Abelian 
functions of time. Integrable tops, as well as almost all classical examples 
from the XlXth century Analytical Mechanics find their place here. 

3. Double loop algebras and their central extensions. This class of algebras 
accounts for integrable PDE's admitting the zero curvature representation 
(such as the Nonlinear Schroedinger equation, the Sine-Gordon equation and 
many others). 

4. Algebras of pseudodifferential operators. The KdV equation and its higher 
analogs come from this example, although it is more practical to derive them 

from double loop algebras."'" 



""^A further generalization is possible: we can add one more spatial variable and consider the 
loop algebra based on the algebra of pseudodifferential operators; this yields the so called KP 
equation for functions of three variables. 
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5. Algebra of vector fields on the line. This algebra, or rather its central exten- 
sion (the Virasoro algebra) and the associated loop algebra again are related 
to the KdV equation. 

As wc sec from this list, the choice of the Lie algebra determines not only the 
possible kinematics (i.e., the structure of the phase space) of the dynamical systems 
which admit a realization based on this algebra, but also the functional class of the 
possible solutions. In all cases, it is very important to examine central extensions 
of the algebras in question (if any), as well as their non-trivial automorphisms: 
they usually lead to new examples. Non- trivial central extensions exist in cases 
2, 3, 4, 5, each one leading to a non-trivial theory. As for outer automorphisms 
of these algebras, they serve, for instance, to define twisted loop algebras with 
important range of applications (see Section 7.1) and are also used in the study of 
nonlinear finite-difference equations. 

It is virtually impossible to provide these lectures with a comprehensive list of 
references. The point of view adopted here is certainly rather subjective. I have 
included a few references to the old original papers as well as to several reviews. 



2. A few Preliminaries: Poisson Brackets, Coadjoint Orbits, 
etc. 



The first question that is prior to the study of the dynamics of a mechanical 
system is that of its kinematics, i.e. of the structure of its phase space. Typically, 
the phase space of an individual dynamical system should be a symplectic manifold. 
However, an important conclusion which may be drawn from the practical study 
of numerous examples is that integrable systems associated with auxiliary linear 
problems always arise in families. The appropriate geometrical setting for the 
study of such families is provided by the theory of Poisson manifolds. Let me 
recall that a Poisson bracket on a smooth manifold M is the structure of a Lie 
algebra on the space C°°(A'/); moreover, the Poisson bracket satisfies the Leibniz 
rule, i.e., it is a derivation with resect to each argument. In local coordinates, a 
Poisson bracket is written as 

{<,,V}(x) = ^7r,,(s)gg, 

where iTij is an antisymmetric tensor (Poisson tensor) satisfying a quadratic dif- 
ferential constraint which assures the Jacobi identity. When the manifold M is 
symplectic, i.e., admits a nondegenerate closed 2-form lu = '^uj^^dxi Adxj, the as- 
sociated Poisson tensor TTy = {uj'-') ^ ■ Reciprocally, whenever the Poisson tensor 
is nondegenerate, its inverse is a symplectic form. In general, a Poisson manifold 
is not symplectic; the fundamental theorem which goes back to Lie asserts that it 
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always admits a stratification whose strata are already symplectic manifolds {the 
symplectic leaves of our Poisson manifold) . The geometrical meaning of this de- 
composition is very simple. Any H € C°° (M) defines on M a Hamiltonian vector 
field which acts onip G C°°{M) via 

XhV = {H, if} ; 

for a given point x G M the tangent vectors Xh (x) span a linear subspacc in 
the tangent space T^M; this is precisely the tangent space to the symplectic leave 
passing through x.^ By construction, Hamiltonian vector fields are tangent to 
symplectic leaves, and hence the Hamiltonian flows are preserving each leaf sepa- 
rately. A closely related property of general Poisson manifolds is the existence of 
Casimir functions. By definition, a function H G C°° (M) is called a Casimir func- 
tion if it lies in the center of the Poisson bracket; equivalently, Casimir functions 
define trivial Hamiltonian equations on M. Restrictions of Casimir functions to 
symplectic leaves in M are constants; reciprocally, the common level surfaces of 
Casimir functions define a stratiflcation of M; typically, this stratification is more 
coarse than the stratification into symplectic leaves (i.e., symplectic leaves are not 
completely separated by the values of the Casimirs), but in many applications the 
knowledge of Casimir functions yields a sufficiently accurate description. 

A very typical example of a Poisson manifold is the dual space of a Lie algebra. 
I shall briefly recall the corresponding construction, since it proved to be very 
important for the study of integrable systems. Let g be a Lie algebra, g* its dual 
space, P {2*) the space of polynomial functions on g* . By the Leibniz rule, a 
Poisson bracket on P (g*) is completely determined by its value on the subspace 
of linear functions g C P (fl*); for X, F G let us set simply 

(2.1) {X,Y}{L) = {L,[X,Y]), Leg*. 

The Jacobi identity for the bracket (2.1) follows from that for the Lie bracket in 
g; since P (g*) is dense in C°° (g*), it canonically extends to all smooth functions. 
Explicitly, we have 

(2.2) {^1, ^2} {L) = {L, [#i {L) , d^2 {L)]) . 

(Note that dipt (L) e (g*)* — g, and hence the Lie bracket is well defined). The 
bracket (2.2) is usually called the Lie Poisson bracket. Its properties are closely 
related to the distinguished representation of the associated Lie group, the coad- 
joint representation. Let G be a Lie group with Lie algebra g; exp : g — > G the 

^The integrability condition which assures the local existence of the submanifold with the 
given tangent distribution immediately follows from the Jacobi identity; the subtle part of the 
proof consists in checking that the possible jumps of the rank of the Poisson tensor tt do not lead to 

singularities of the leaves; one can check that the Lie derivative of tt along any Hamiltonian vector 
field is zero, and hence its rank is constant along the leaves (but may jump in the transversal 
direction) . 
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exponential map. The adjoint and coadjoint representations of G acting in g and 
0*, respectively, are defined by 



Adg-X = (^j^ g-exptX -g ^, X €Q, 
{Ad*g-L,X) = {L,Adg-^-X), Xgs, Leg*. 



Set 



adX ■ Y = (4-] Ad exptX -Y, ad*X-L = ( 4 ) Ad* expiX • L. 

Clearly, one has adX-Y = [X,Y] , ad* X = — (adX)* . The following fundamental 
theorem again goes back to Lie; it was rediscovered by Kirillov and Kostant in 
1960's: 



Theorem 2.1. (i) Symplectic leaves of the Lie-Poisson bracket coincide with G- 
orbits in g* (coadjoint orbits), (ii) Casimir functions of the Lie-Poisson bracket 
are precisely the coadjoint invariant functions on q* . 
It is very easy to verify a somewhat weaker property. 



Proposition 2.2. Let (f G C°°{q*) be an arbitrary function; the Hamiltonian 
equation of motion defined by ip with respect to the Lie-Poisson bracket may be 
written in the following form: 

(2.3) ^ = ~ad*d(p{L)-L, Lgq*; 

in other words, the velocity vector, associated with any Hamiltonian equation on 
Q* is automaMcally tangent to the coadjoint orbit passing through L. 

In the context of integrable systems coadjoint orbits arc of particular importance: 
in many applications, the phase spaces of integrable systems are coadjoint orbits 
for some appropriate Lie group^. Classification of coadjoint orbits for particular 
Lie groups is a good exercise (which may be quite involved depending on the 
nature of the Lie group); let us just quote a few examples which will be useful in 
the sequel. 

^It's probably worth making some precisions: when we deal with concrete equations, coadjoint 
orbits are almost allways a good starting point, but it may be quite useful to enrich our tools. 
In some cases, an orbit is too big for our purpose and it is possible to cancel out some degrees 
of freedom by passing to the quotient space over some manifest symmetry group. On the other 
hand, in some cases, an orbit is too small and it's more practical to use a bigger phase space 
which is mapped onto the orbit in a way which is compatible with its Poisson structure. Finally, 
there arc classes of examples when the context of Lie algebras appears to be too restrictive and 
we have to deal with nonlinear Poisson brackets from the very beginning. We shall comment on 
these examples later on (see Section 9). 
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Example 2.3. Let g = Ql{n) be the full matrix algebra; its dual space g* may be 
canonically identified with Q by means of the invariant inner product^ 

(2.4) {X,Y)=trXY. 

Thus the adjoint and the coadjoint representations of the corresponding Lie group 
G = GL{n) are identical; we have Ad*g-L = gLg~^; the coadjoint orbits consist of 
conjugate (isospcctral) matrices; their classification is given by the Jordan normal 
form. Casimir functions are spectral invariants of matrices; their level surfaces 
consist of a finite number of coadjoint orbits. 

Example 2.4. Let b+ C be the subalgebra of upper triangular matrices; the 
pairing (2.4) allows to identify its dual with the space b_ of lower triangular 

matrices. The coadjoint representation of the corresponding Lie group of 
upper triangular invertible matrices is given by 

Ad*b-F = P_{b-F-b-^), b G B+, Feb-, 

where P_ : g ^ b_ is the projection operator which replaces by zeros all matrix 
coefficients above the principal diagonal. 

In this example the adjoint and the coadjoint representations are inequivalent. 

After this brief discussion of Poisson geometry and coadjoint orbits let us return 
to the study of integrable systems. The use of linear Poisson brackets and of 
coadjoint orbits seems a good guess to get a proper kinematical description of 
our future examples; this suggestion is further supported by the following simple 
observation which specializes proposition 2.2 above. 

Proposition 2.5. Assume that q = flt(n) is identified with its dual space and 
equipped with the Lie-Poisson bracket. For any G C°° (g) the Hamiltonian equa- 
tion of motion is written in the form 

(2.5) ^ = -[d^{L),L]; 

hence all Hamiltonian flows on g preserve spectral invariants of matrices. 

Equation (2.5) looks exciting: one is tempted to compare it with the famous 
Lax equations. A closer look on the picture reveals, however, that proposition 2.5 
leads to a deception. Indeed, spectral invariants of matrices are Casimir functions 
for the Lie-Poisson bracket; their conservation is a trivial fact which has nothing 
to do with integrability of equation (2.5). There is very little chance that this 
equation with arbitrary Hamiltonian ip will be completely integrable; on the other 
hand, tlic sptx tral iiuariauls themselves which seem to be natural candidates to 

^Throughout these lectures inner product means a nondcgcncrate symmetric bilinear form 
(C- bilinear in the case of complex algebras); over the reals we do not impose any positivity 
condition. 
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provide integrable systems, generate trivial flows, in view of the following simple 
fact: 

Proposition 2.6. For any Lie algebra q a function e C°°(fl*) is a Casimir 
function for the Lie-Poisson bracket on q* if and only if 



for any L £ q* . (Note that dip (L) g (g*)* — Q; when q and q* are identified, this 
relation is reduced to [dip (L) , L] = 0.) 

Despite this initial setback, the original idea to use Lie-Poisson brackets and 
coadjoint orbits can be saved. However, instead of the initial Lie-Poisson bracket 
which provides the set of spectral invariants but does not yield any nontrivial 
dynamics associated with them, we must find a different one. It's at this point 
that the classical r-matrix is brought into play. 



3. Classical r- matrices and LcLx Equations 

Let g be a Lie algebra We shall say that r G End g is a classical r-matrix if the 
bracket 



is a Lie bracket, i.e. if it satisfies the Jacobi identity. The skew symmetry of (3.1) 
is obvious for any r. Wc denote the Lie algebra with the bracket (3.1) by fl^ and 
say that {g, Qr) is a double Lie algebra. 

If is a double Lie algebra, there are two diflferent Poisson brackets in the space 
g*, namely, the Lie-Poisson brackets of g and 0r- The latter bracket will be referred 
to as the r-bracket, for short. 

A class of double Lie algebras which is of great importance for applications is 
constructed as follows. Assume that there is a vector space decomposition of g 
into a direct sum of two Lie subalgebras, g = 0+ + Let P± be projection 
operators onto g± parallel to the complementary subalgebra; set 



ad* d(p (L) ■ L = 



(3.1) 



[X,Yl = -{[rX,Y] + [X,rY]) 



(3.2) 



r = P+-P_. 



In this case, bracket (3.1) is given by 



(3.3) 



[X,Y]r = [X+,Y+[ 



] - [X-Y.] 



where X± = P±X,Y± = P±Y. In other words, the bracket (3.1) is the difii'erence 
of Lie brackets in 0+ and 0_. The Jacobi identity for g^ is obvious from (3.3). 
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As discussed in Section 6, in typical applications the Lie algebra g is a loop 
algebra, i.e., an algebra of matrix- valued functions on the circle, and q± its subal- 
gebra consisting of functions which are analytic inside (resp., outside) the circle. 
In that case, the classical r-matrix (3.1) is precisely the Hilbert transform. Of 
course, general classical r-rnatriccs need not have this simple form (although (3.1) 
is by far the most important example of all). We shall discuss the general theory 
of r- matrices a little later; let us first state the key theorem which motivates the 
definition. 

3.1. Involutivity Theorem 

Let I{g*) be the ring of Casimir functions on g* (with respect to the original Lie- 
Poisson bracket); equivalently, C C°°{g*) is the set of coadjoint invariants. 

Theorem 3.1. (i) Functions in I{q*) are in involution with respect to the r- 
bracket on q* . (ii) The equations of motion induced by h & I{q*) with respect to 
the r-bracket have the form 



If Q admits a nondegenerate invariant bilinear form, so that ad* ~ adg, equations 
(2.4) have the Lax form 



Proof, (i) Let /ii,/i2 G /(fl*); set dht {L) = Xi for brevity. By definition. 



(3.4) 



— = -ad^M -L, M = r{dh{L)) 



{huh2}, {L) = {L,[Xi,X2l) 

= ^{L,[rX,,X2] + [X,,rX2]) 

= ^ {ad;X2 ■ L, rXi) - i (ad*Xi • L, rXs) = 0, 



since, by proposition 2.6, ad*X2 ■ L = ad*Xi ■ L = 0. (ii) We have 



— - = —ad* dh (L) ■ L; 



(3.1) implies that 



(3.5) 
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since h G /(fl*), the second term in (3.5) vanishes. 



□ 



Remark 3.2. The matrix M in (3.4) is not uniquely defined: one can always add 
to it something which commutes with L. Here is a useful option: set 



equation (3.4) holds with any of these two operators. Below, we shall see that this 
form of M-operator appears naturally from the global formula for the solutions. 

Theorem 3.1 has a transparent geometrical meaning: it shows that the trajecto- 
ries of the dynamical systems with Hamiltonians h £ I{q*) he in the intersection 
of two families of orbits in g*, the coadjoint orbits of g and g^- Indeed, the coad- 
joint orbits of Qr are preserved by all Hamiltonian flows in g^- On the other hand, 
because of (2.4), the flow is always tangent to the g-orbits in g*. 

In many applications the intersections of orbits are precisely the "Liouville tori" 
for our dynamical systems. 

3.2. Factorization Theorem 

The scheme outlined so far incorporates only two of the three main features 
of the inverse scattering method: the Poisson brackets and the Lax form of the 
equations of motion. As it happens, the most important feature of this method, the 
reduction of the equations of motion to the Riemann problem, is already implicit 
in our scheme. An abstract version of the Riemann problem is provided by the 
factorization problem in Lie groups. 

We shall state a factorization theorem, which is the global version of Theorem 
3.1, for the simplest r-matrices of the form (3.2). Let G be a connected Lie group 
with Lie algebra g, and let G± be its subgroups corresponding to g±. 

Theorem 3.3. Let h E I (g*) ,X = dh (L). Let g±{t) be the smooth curves in G± 
which solve the factorization problem 



Then the integral curve L{t) of equation (2.4) with L{0) = L, is given by any of 
the two formulae, 



(3.6) 




(3.7) 



expiX = g+{t) ■ g-{t)-\ g+{0) = e. 



(3.8) 



L{t) = Ad*Gg+{t)-' ■ L = Ad*Gg-{tr' ■ L. 



Proof. Differentiating (3.8) with respect to t we get 



dL 
dt 



ad* {g+^g+) ■ L = -ad* {gJ^g-) ■ L. 
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We shall check that g±^g± = M±, where M± are the M-operators from (3.6). Due 
to our special choice of r we have 

M± it) = ±P±X {t) , 

where X (t) = dh {L (t)). The Ad* G-invariancc of h implies that 

x{t) = AdGg±{t)-'-x. 

Exercise 3.4. Check this formula in matrix case ("gradients of invariant functions 

are covariant"). 

Writing (3.7) in the form g+(t)exptX — g-(t) and differentiating with respect 
to t, we get 

g+^g+ + AdGg-{t)~^ ■ X = g-gZ^ 
Since g±^g± & Q±, this implies g±.^g± = ±P±X {t), as desired. □ 

Note that the two possible choices of sign in (3.8) are equivalent precisely because 
X{t) belongs to the centralizer of L{t), i.e., ad*X{t) ■ L(t) = (in fact, this is the 
characteristic property of Casimir functions). 

By the implicit fimction theorem, the factorization (3.7) exists for t sufficiently 
small; note that in our proof we need not assume that this factorization exists 
globally for all t. Geometrically, this means that the solution of the Lax equation 
exists as long as the curve exp tX remains in the "big cell" G+ x G_ C G; in other 
words, the flow associated with the Lax equation is not necessarily complete. 
One can show in typical examples that the curve intersects "complementary" cells 
of positive codimension transversally and returns back to the big cell: for the 
exceptional values of t the solution "escapes to infinity", i.e., displays a pole in t. 

3.3. Factorization Theorem and Hamiltonian Reduction 

A more geometric proof of Theorem 3.3 is based on Hamiltonian reduction. 
Recall that the Hamiltonian reduction applies to Hamiltonian dynamical systems 
with high degree of symmetry; it allows to exclude certain redundant degrees of 
freedom. Classically, the use of reduction is to simplify multidimensional systems 
getting quotient, or reduced, systems of lower dimension. However, as pointed out 
in [9] , one can reverse this reasoning and use Hamiltonian reduction in the opposite 
direction, starting with a simple multidimensional system with high symmetry 
("free dynamics") and getting a complicated reduced system as an output. In 
order to apply this idea, one has to answer the following questions: 

1. Find a "big" phase space and suggest the "free dynamics" which will yield 
Lax equations as the quotient system. 

2. Make sure there is an expected high degree of symmetry for the free system. 
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3. Perform the reduction. 

Although the proof based on this approach is much longer than the elementary 
computation presented above, it is more transparent and explains the origin of the 
result. A simple candidate for the big phase space is the cotangent bundle T*G 
equipped with the canonical symplectic structure. Let us first of all describe the 
"free dynamics" on T*G. The group G acts on itself by left and right translations; 
these actions naturally lift to T*G', both actions arc Hamiltonian with respect to 
the canonical symplectic structure. Let us identify T*G with G x g* by means of 
left translations. 

Exercise 3.5. In left trivialization the action of G on T*G ~ G x by left (right) 
translations is given by 

(3.9) \{g) : (x, L) ^ (gx, L) , 

p{g) : {x,L)^ {xg'\Ad*g-L) . 

(In left trivialization the action of G by right translations induces a nontrivial 
action in the fiber g*; it is easy to check that it is precisely the coadjoint ac- 
tion.) Left-invariant functions on T*G are identified with fimctions on g*. Since 
the canonical Poisson bracl?;et of left- invariant functions is also left-invariant, this 
induces a Poisson structure on g*; it is easy to check that it coincides with the 
Lie-Poisson bracket. Casimir functions on g* canonically lift to smooth functions 
on T*G which are G-bimvanant. For h G /(g*) let us denote by ft € G~ {T*G) 
the corresponding biinvariant Hamiltonian on T*G. 

Lemma 3.6. The Hamiltonian flow on T*G defined by h is given (in left trivial- 
ization) by 

(3.10) Ft : {x, L)\ — >{x- exptdh{L), L), xgG, Leg*. 

In other words, integral curves of h project to left translates of one-param,eter 
subgroups in G; the choice of h determines the (constant) velocity vector dh{L) 
which depends only on the initial data. 

Since the "free Hamiltonian" h is G-biinvariant, the flow Ft admits reduction 
with respect to any subgroup U C G x G. There is, at this stage, a very large 
freedom in the choice of such a subgroup which all lead to different but meaningful 
quotient systems. The particular choice which is imposed by the choice tlic r- 
matrix (3.2) is U = G+ x G_. By (3.9), with our choice of trivialization, the 
action of G+ x G_ on T*G ~ G x g* is given by 

(3.11) (.g+,.g_) : {x,L) ^ {g+xgZ\ Ad},g_ ■ L) . 

We now turn to the reduction procedure. In textbooks the reduction is usually 
described in a rather complicated way (which involves the moment map and a good 
deal of symplectic geometry) (see, e.g., [7]). Here is an elementary substitute. Let 
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us ask what is bad about the naive suggestion: consider a group action GxM —>■ M 
on a symplectic manifold and take the projection -k : M M/G onto the quotient 
space? The answer is: the quotient space is no longer symplectic, since symplectic 
forms transform by pullback, and there is no natural symplectic form on M/G (it 
in not even in general even dimensional!) But on the other hand, the quotient 
space does carry a Poisson bracket (which transforms by push-forward!). The 
difficult part of reduction consists in the description of the particular symplectic 
leaves of this quotient Poisson bracket; it's at this stage that one needs the moment 
map and all other machinery. If we do not want a too detailed description, or if 
we manage to guess the symplectic leaves in some other way, everything becomes 
simple! In the present case, we can display a map tt which is constant on the orbits 
of G-i- X G- in M and hence its image yields a model of the quotient space; the 
Poisson structure on this quotient is easy to recognize. 
For X G G we denote by x± the solution of the factorization problem 



(3.12) 



: x+ ■ x_^, x+ € G+, X- e G- 



L is constant on 



Lemma 3.7. (i) The map tt : T*G — > g* : [x, L) i — > Ad}.: 
the orbits of G+ x G- in T*G. (ii) If G is globally dijfeomorphic to G+ x G-, i.e., 
the factorization problem (3.12) is always solvable, it is a global cross-section of 
this action. (Hi) The induced Poisson structure on T*G/G+ X G- ~ Q* coincides 
with the Lie-Poisson bracket for — 0+ © 0-- 

The check of (i) and (ii) is immediate; the proof of (iii) requires a little knowledge 
of symplectic geometry (or else a three-line computation); we shall not present it 
here (see [8]). 

Lemma 3.8. The flow (3.10) factorizes over q*; in other words, there exists a 
natural flow Ft : Q* — > Q* (called the quotient flow) which makes the following 
diagram 



T*G 



T*G 



commutative; the quotient flow Ft : g* — > is given by 



(3.13) 



Ft:L^Ad*G9- {ty'-L 
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where, as in (3.7), g+ (t) 5_ (t) = exptdh{L). 

The flow Fi is precisely the result of the reduction procedure. Note that (3.13) 
involves only g_; this is due to our choice of trivialization of T*G; trivialization 
by right translations yields the equivalent formula for the quotient flow Ft : L i — > 
Ad*ag+ {t)-'-L. 

In general, of course, G need not be diffeomorphic to G+ x (?_; still, Q* provides 
a model for a "big cell" in the quotient space T*G/G+ x G„; one can show that 
under very mild restrictions the action (3.11) is proper and hence the quotient 
space is a well-defined manifold; the quotient flow induced on this manifold is also 
well defined and may be regarded as the natural completion of the incomplete fiow 
associated with Lax equations. 

The choice of the subgroup G+ x G_ C G x G as the symmetry group for the 
"free system" may seem arbitrary; indeed, there are many other possible choices 
leading to meaningful examples. (Among the dynamical systems which may be 
obtained in this way, there are, e.g., the Calogero-Moser systems, cf. [9].) The 
key property which characterizes our special situation is the simple description of 
the quotient space: as we see from Lemma 3.8, in this case the quotient space is 
simply the dual of a Lie algebra with its Lie-Poisson bracket, and its symplectic 
leaves (which are the symplectic quotients of T*G) are the coadjoint orbits of G^; 
for other choices of the symmetry group the description of the quotient space will 
be more complicated and it may no longer be a homogeneous space. 

4. Classical Yang-Baxter Identity 

We have already mentioned that the most interesting r-matrices are associated 
with decompositions of the Lie algebra into complementary Lie subalgebras. How- 
ever, it is worth examining the general conditions on r which follow from the Jacobi 
identity for the r-bracket. These conditions are known as the classical Yang-Baxter 
equations; they were first derived as a semi-classical approximation to the quantum 
Yang-Baxter equations which arise in the study of quantum completely integrable 
systems. 

The restrictions on r which follow from the Jacobi identity are quite easy to 
establish. For r € End g set 

(4.1) Br {X, Y) = [rX, rY] - r{[rX, Y] + [X, rY]). 

Proposition 4.1. The r-bracket (3.1) satisfies the Jacobi identity if and only if, 
for any X,Y,Z G g, 

(4.2) [Br{X,Y),Z] + [Br{Y,Z),X] + [Br{Z,X),Y] = 0. 

The proof is straightforward: just substitute (3.1) into the Jacobi identity and 
regroup the terms. The necessary and suSicient condition (4.2) is usually replaced 
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by sufficient conditions which are bilinear rather than trilinear. The simplest 
sufficient condition is the so-called classical Yang-Baxter equation (GYBE) 

(4.3) Br{X,Y)=0. 

Another important sufficient condition is the modified classical Yang-Baxter equa- 
tion (mCYBE) 

(4.4) Br (X, Y) = -c[X,Y], c = const. 

By rescaling, we may always assume that c = 1. Note that the r-matrices (3.2) 
satisfy mCYBE with c = 1. 

The reason for the study of classical r-matrices satisfying the modified classical 
Yang-Baxter identity is that although they do not in general have the simple form 
(3.2), one can still associate with them a factorization problem. By contrast, the 
ordinary classical Yang-Baxter identity (4.3) represents a degenerate case and does 
not lead to a factorization problem. 

Let us briefly describe the corresponding construction. Given an r-matrix which 
satisfies mGYBE, set 

(4.5) r± = i (r ± Id) . 

Proposition 4.2. For each X,Y £ q 

[r±X,R±Y]=r±{[X,Y]r), 

i.e., r± : Qfj g are Lie algebra homomorphisms. 

Set g± — Imr±. Clearly, g± is a Lie subalgebra of g. If r has the form (3.2), 
then r± = ±P± and the subalgebras q± = P±{2) are complementary. In general 
case this is no longer true. However, this difficulty may be resolved by passing to 
the double of q. By definition, 

0=000 

is the direct sum of two copies of 0. 

Proposition 4.3. The mapping ir : 2r ^ g (B Q : X i — > (r+X, r_X) is a Lie 
algebra embedding, and each Y € g has a unique decomposition, Y = Y^ —Y-, 

where (y+ , y_ ) G Im i^ . 

The last assertion follows immediately from the obvious identity r+ — r_ = Id. 

Now, let G, Gr be (local) Lie groups which correspond to 0, 0^. The homomor- 
phisms r± give rise to the Lie group homomorphisms which we denote by the same 
letters. Put G± = r±(G_R). The composition of maps 

V : Gr — >■ G X G : X i — »■ {r+x, r-x) 

is a Lie group embedding. Gonsider the map 

m : G X G — > G : {u, v) i — > uv~^. 
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Then f = moi^ : Gr 



G is a local homeomorphism, and therefore an arbitrary 



element y G G which is sufficiently close to unity admits a unique representation 

y = y+yZ^ with {y+,y-) e Imir. 

The proof of Theorem 3.3 extends to the present setting with only minor changes. 

It is probably worth giving examples of r-matrices satisfying (4.4) which are 
not of the form (3.2). Let g = fll(n) be the matrix algebra; let us consider its 
decomposition q = n+-i-f)-t-n_ into direct sum of upper triangular, diagonal, and 
lower triangular matrices. Let P±,Po be the corresponding projection operators. 
Let ro be the partially defined linear operator on q with domain Dom(ro) = n+-i-n_ 
given by 



We want to extend ro to the entire linear space g in such a way that it satisfies 
the modified classical Yang-Baxter identity. Let us first drop the latter condition 
and consider all linear operators r D ro, Dom(r) = g. Let us set again r± = 
^ {r ± Id) ,ir = r+ ® r_ and consider the subspace Qr = Imi^. It is easy to 
see that this subspace is transversal to the diagonal subalgebra Qd = {{X,X) e 
0®£l}i conversely, all extensions r D ro are in bijective correspondence with linear 
subspaces flr C g © which contain n+ ® n_ and are transversal to the diagonal. 
All such subspaces may be parametrized in the following way: 



where X± = ±P±X, Xq = PqX and G End [) is a linear operator; equivalently, 
the extensions r D ro are described by the von Neumann formulae, 



The transversality condition is equivalent to the non-degeneracy of Id — 0. Since 
f) C 5 is an abelian subalgebra which normalizes both n+ and n+, it is easy to 
see that flr C ® is a Lie subalgebra (and not merely a linear subspace) for 
any 9, and hence all re D ro satisfy mCYBE. In more general examples, when the 
relevant subalgebra is not abelian, there are additional algebraic constraints on 6 
which make the construction partially rigid. 

The moral of the story: One can construct new classical r-matrices from the 
standard ones of the form (3.2) by assuming that the simple formula (3.2) ap- 
plies not globally, but only on some subspace of positive codimension and then 
using the extension theory of linear operators. The same trick works for loop al- 
gebras; in that case, the relevant codimensions are finite. Let us set, for example, 
25 = Lg, Q = 0[(n); by definition, © is the algebra of Laurent polynomials with 
matrix coefficients. Set = {X{z) = Xq + Xiz + X2Z^ + ...]Xo e n+}, 0T_ = 
{X{X) =Xo + Xiz-'^+X2Z-'^ + ...;Xo € n_}. Then (5 = 01++f)-i-9l_, where f) is 
the subalgebra of constant diagonal matrices. Hence we may apply the previous 




0^ = {{X+ + Xo, X- + 0Xo) ■,X€g} 
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construction without any modification. A classification theorem, due to Belavin 
and Drinfeld, assures that all r-matrices on semisimple Lie algebras and their loop 
algebras satisfying mCYBE and certain additional conditions arise in a similar 
way; of course, the most difficult part of the construction is the case when f) is 
replaced by a non-abelian subalgebra. (Cf. [10], [11].) 

5. A Finite-dimensional Example 

Let again g = gl{n) be the full matrix algebra. There are several natural decom- 
positions of Q into direct sum of complementary subalgebras, e.g., 

(5.1) = b+ + n_, 

where n_ is the subalgebra of lower triangular nilpotent matrices and b+ = f) + n+ 
the complementary subalgebra of upper triangular matrices. Another decomposi- 
tion is 

(5.2) fl = b+ + fi, 

where 6 = so{n) is the Lie algebra of skew symmetric matrices. We can associate 
two standard classical r-matrices with these decompositions: 

^Gauss — -^n- ; 

^Iwasawa — 

where Pb+,Pn-,Pb+,Pt are the respective projection operators (mind that of 
course Pb^ 7^ Pb+)- The Lie groups Grc^^^^ and Griwasawa are isomorphic to 
X A''_ and i?+ x K, respectively. The associated factorization problems in 
the general linear group G = GL{n) are the Gauss decomposition of matrices, 

(5.3) g = bn~^,b€ B+,nG N-, 
in the former case, and the Iwasawa decomposition 

(5.4) g = bk-\beB+,keK, 

in the latter one. (Here N- c GL{n) denotes the subgroup of lower triangular 
unipotcnt matrices, C GL{n) the subgroup of upper triangular matrices, and 
K = SO{n) C GL(n) the subgroup of orthogonal matrices.) Note that in the 
Iwasawa case the product map _B+ x K ^ G : {b,k) ^ bk~^ is a bijection onto G 
and hence the factorization problem is always solvable. For the Gauss decompo- 
sition the image of B+ x iV_ G : {b, n) i — > bn~^ is an open dense subset in G, 
and hence the factorization problem is solvable for almost all (though not for all) 
initial data. 



18 



Semenov-Tian-Shansky 



The dual space q* is canonically identified with g by means of the invariant inner 
product 

{X,Y) = trXY; 

the decompositions (5.1), (5.2) give rise to the biorthogonal decompositions 

fl* = b| + n:^ = b| + e^, 
which provide the models for dual spaces of the subalgebras, 

~ n:): = b-,n*_ ~ b:}: = n+ 
in the case of the Gauss decomposition and 

b;:^«^ = p,r ~bi = n+ 

in the case of the Iwasawa decomposition (here p C Mat{n) denotes the subspace 
of symmetric matrices).^ Let us focus on coadjoint orbits of B_^_. Denote by 
dp C Mat{n),p e Z, the set of all matrices supported on p-th diagonal, 

j)p = {X G Mat{n);Xij = for j - i ^ p} 

Exercise 5.1.(i) Let us model the dual space b^j. on b_; for all g > the subspaces 
b'L = ®p^o^-p C ^- invariant with respect to the coadjoint action of B+. (ii) 
Similarly, if b^ is modelled on p, the subspaces p' = ®p=o i^p + ^-p) ^ P 
invariant. 

Orbits in the subspace b^ are all trivial (i.e., each point in this subspace is stable 

with respect to the coadjoint action and hence is a separate orbit); the subspace 
bL contains orbits of maximal dimension 2n — 2 (over C there is just one such 
orbit, over the reals there is a finite number of them); the typical example is the 
orbit Of which contains the matrix /el)_iCb_, 

/ ••• \ 

: "■• 

Vo •■• 1 0/ 



Exercise 5.2. (i) Of consists of all matrices of the form 
/ pi • • • \ 

6i ■•• ■•• ; 

: '■• Pn-l 



/ = 



,&i?^0,^ft = 0. 



Pn } 



^Mind that the model for the dual of a subalgcbra depends not only on the subalgcbra itself 
but also on the choice of its complement in the big Lie algebra (and eventually also on the choice 
of the inner product whenever it is not unique). 
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(Over the reals, in addition, signbi = +1.®) (ii) Lie-Poisson brackets of coordinate 
functions Pi, bj on Of are given by 

{PiiPo} = {bi,bj} = 0, {pi,bi} = - {pi+i,bi} = bi. 

(iii) If we set bi = exp (qi — qi+i), the coordinates qi have canonical Poisson brack- 
ets with momenta, {pi,qj} = Sij. 

Thus, as a symplcctic manifold. Of is isomorphic to the standard phase space 
-g2n-2^ In the Iwasawa model, the same orbit is realized by symmetric matrices. 



(5.5) 



/ Pi bi 
bi ■■• 

V ■•• 





Pn-1 bn-1 
bn-1 Pn 



\ 



Recall that our main theorem associates dynamical systems to coadjoint orbits in 
fl^; since Qr^^^, and Qr 



^Iwasawa 
fl'"Gauss 



are direct sums. 



respectively. Coadjoint orbits of this bigger algebra are Cartesian products of the 
coadjoint orbits of the factors. In the Iwasawa case the simplest meaningful choice 
is to take the zero orbit of K in ~ n+; with this choice, (5.5) becomes the 
Lax matrix for the associated dynamical systems. The Hamiltonians are spectral 
invariants of L; taking, for instance, 

2 

and expressing it in terms of the canonical variables Pi,qi, we get 
^ = J I^P' + I] exp 2 {qi - qi+i) , 



i.e., the Hamiltonian of the open Toda lattice. In the Gauss case, it is also possible 
to take the zero orbit of the complementary subalgebra n_ ; with this choice the 
Lax matrix I will be lower triangular and hence its spectral invariants will depend 
only on momenta Pi and the corresponding Hamiltonians will be trivial. Luckily, 
in this case there is another option: we may take any one-point orbit {e} of N- 
in n1 — n+; the constant matrix e is simply added to the Lax matrix. This 
pro( 0(hii o (loos not add an\ new degrees of freedom to our system, but it modifies 



^It is easy to see that in the real case the signs of all matrix coefficients below the principal 
diagonal are preserved by the coadjoint action; hence there are exactly 2"~^ open orbits in our 
subspace. 
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the embedding of the "httle orbit" Of into the big algebra and hence the spectral 
invariants of the Lax matrix. Specifically, set 



/ 1 



\ 



\ 



1 

/ 



/ + e. 



Remark 5.3. The reader familiar with semisimple Lie algebras will notice that the 

Jordan matrices e, / arc principal nilpotent elements of the general linear algebra, 
which suggests the way to generalize the above construction to other algebras. 



Exercise 5.4. (i) {e} C 5i C n+ is a one-point coadjoint orbit of 

(ii) In canonical coordinates, the Hamiltonian h = ^tr Lq^^^^ is given by 

i i 

Note that the Hamiltonians H and h (which are defined on the same manifold 
Of) arc different (although, in this particular case, they happen to be related by a 
simple canonical change of variables); so are the associated factorization problems 
which solve the Toda equations. 

Other spectral invariants, e.g., hp = L^^^^^^ P = 2,3,..., form a system of 
integrals of motion in involution. Obviously, the number of independent integrals 
does not exceed n — 1 (mind that hi reduces to a constant on Of). 



Exercise 5.5. (i) Write down explicitly the two factorization problems, (ii) Find 

the explicit relation between them, (iii) Show that over the reals the group element 
exptLcaussiP, q) lies in the "big cell" for all (p, q) G M^""^ and for all t €R, and 
hence the Gauss factorization is always possible, (iv) Check that the integrals 
h2, ■■■,hn remain fimctionally independent after restriction to Of. 

In both cases, the entries of the factors are rational functions of exp t (with coef- 
ficients depending on the initial data, i.e., on p, q)\ thus the functional dependence 
of the solutions on the time variable is fairly simple. It is easy to see that this sim- 
ple behaviour is characteristic for all Lax equations associated with factorization 
problems in finite-dimensional Lie groups. 



Remark 5.6. The subspaces £)p C g satisfy [Op, Dq] C fp+q, i.e., they form a grad- 
ing of the matrix algebra. As a consequence, the subspaces ®p>n5p, n = 0, 1, 2, 
are Lie subalgebras and form a decreasing filtration of b+; by duality, the subspaces 
= ©p<_„Z)_p C b_ ~ b!j_, n = 0,-1,-2,..., form an increasing filtration of 
b\ by Ad -invariant subspaces. The orbit Of is the "biggest" orbit in the sub- 
space Oi; its choice is quite natural, since dimO/ = 2n — 2 is twice the number 
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of independent Hamiltonians, which is precisely the amount needed for complete 
integrability. Dynamical systems associated with coadjoint orbits of higher dimen- 
sion (which are abundant) have not got enough "obvious" integrals of motion to 
assure their Liouville integrability; on the other hand, all these systems are explic- 
itly solvable by means of the factorization problem. This queer situation is due to 
the resonance behaviour of these systems: their trajectories span a submanifold in 
the phase space with codimension higher than in the generic case. 

An important general conclusion to be held from this example is the role of 
grading: it provides a natural decomposition into subalgebras, as well as plenty of 
invariant subspaces for the coadjoint action. 



6. Loop Algebras and the Riemann Problem 

As already mentioned, loop algebras provide a natural environment for the study 
of numerous finite-dimensional systems. In this Section we shall briefly outline the 
corresponding constructions. 

Let be a semisimple Lie algebra with invariant inner product (g = sl(n) with 
the inner product (X, Y) = ti XY is a good example; in the sequel we shall mainly 
deal with this standard matrix case). Its loop algebra, Lg is the Lie algebra of 
Laurent polynomials with coefficients in g, Lq — 2 [z, z^^'\ with pointwise commu- 
tator, [Xz'',Yz'^] = [X,y]^"+™, or [x,y] (z) = [x{z),y{z)]. We may regard an 
element x G (5 as a polynomial mapping from the unit circle T into g. An invariant 
inner product on © is given by 

f dz 

(6.1) {x,y) = j^tvx{z)y{z)^^. 

Let Qn = 0-2;"; clearly, [0n,0in] C Qn+m, and hence the decomposition 
Lq = Qn defines a grading (the so called standard grading) of the loop 
algebra. Set 

(6.2) Lq+ = e„>0fln = M , = ®n<0Qn = Z~'^Q [z''^] . 

Proposition 6.1. (i) Lqj^ and Lq_ are graded subalgebras of Lq, and 

Lq = Lq++Lq_. 

(ii) The inner product (6.1) sets Lq into duality with itself; in particular, Lq± is 
set into duality with I/flq:. (Hi) The ring of Casimir functions on Lq* ~ Lq* is 
generated by the functionals 

^n,m [L] = -Res,=otrL(z)"z-"-i. 
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Note that these functionals are smooth in the sense of ordinary calculus of vari- 
ations and grad$„ „i — -^L^~^. 

Let P± be the projection operators onto Lq± parallel to the complementary sub- 
algebra, r = P+ — P_. In analytic terms, Lq consists of trigonometric polynomials 
on the circle, and r is the standard Hilbert transform. The Lie algebra {Lq)^ is 
isomorphic to the direct sum Lqj^ ® Lq_. For X G we set X± = ±P±X. 

Proposition 6.2. The coadjoint action of (Lg)^ on its dual is given by 

ad*X ■ L = P+[X_,L]- P_ [X+, L] ; 

this action leaves invariant the subspace of polynomial loops Lg C (Lg)^. 

Notice that both linear operators L i-^ and L i-^ P_ [X_|_,L] are 

Toeplitz. 

Since Lg is infinite-dimensional, the choice of the associated Lie group becomes 
non-obvious. A reasonable option is to take the group Gw of all Wiener maps 
g -.T ^ G. The Lie algebra Lg may also be replaced by its appropriate completion, 
e.g. the Wiener algebra Lgw- Of course, with this choice the full dual space Lg^,, 
becomes a rather complicated object; the point is that the set of polynomial loops 
Lg c Lg^ is invariant with respect to the coadjoint action of (Lg)^ (though of 
course not with respect to the coadjoint action of Lgw^-)- Let C/^, C Qw be the 
subgroup of Wiener maps (7 : T — > G which are holomorphic in the unit circle, and 
Qw C Gw the subgroup of maps which are holomorphic outside the unit circle 
and satisfy the normalization condition lim^^oo 9 (z) = id. The Lie group which 
corresponds to {Lg)^ may be identified with Gly = Gw x Gw t its coadjoint action 
is given by 

{Ad*h-L) = P+{h^{z)L{z)hZ^ {z)) - P-{h+{z)L{z)hl^ {z))] 
h = {h+,h-) eG^ xGw- 

Exercise 6.3. For g = sl{n) describe all coadjoint orbits of Qw in the subspace 
of Laurent polynomials of the form 

L{z) =l-iz-'^ +lo + hz. 

Remark 6.4. An analyst may feel disappointed with our choice of the 'restricted 
dual' Lg C ifl^ of the Wiener algebra; indeed, the coadjoint orbits which are 
contained in this space are modelled on (matrix) polynomial functions; even gen- 
eral rational functions are not allowed, to say nothing of more interesting classes 
of analytic functions. The reason for this deliberate restriction is very simple: we 
are willing to get dynamical systems which admit a simple parametrization and 
(possibly) some physical interpretation; practical experience shows that examples 
which are physically interesting are usually associated with coadjoint orbits of the 
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lowest possible dimensions. (By contrast, most of the orbits which lie in the 'ex- 
otic' part of the full dual space are infinite-dimensional.) This does not mean of 
course that dynamical systems which are modelled on analytic functions of a more 
complicated nature are totally uninteresting; but once again, 'good' examples are 
associated not with generic coadjoint orbits in the 'very big dual', but rather with 
well-embedded finite-dimensional ones. (Below we shall see how to construct such 
examples using a different choice of the basic Lie algebra.) One more reason to 
single out the finite-dimensional orbits is the possibility to bring into play the 
highly powerful machinery of Algebraic Geometry: wc shall sec below that poly- 
nomial (or, more generally, rational) Lax matrices give rise to algebraic curves of 
finite genus and Lax equations are linearized on their Jacobians. Lax matrices 
associated with infinite-dimensional orbits will lead to curves of infinite genus. 

The specialization of our main theorem to the present situation may be stated 
as follows: 



Theorem 6.5. (i) Invariant junctionals <&„,rn, give rise to Hamiltonian equations 
of motion on Lq c (ig)* with respect to the Lie-Poisson bracket of Lq^; these 
equations may be written in the Lax form, 

(6.3) ^ = [L, M±] , M± = ±P± (grad $„,^) . 

(ii) The integral curve of (6.3) with origin Lq is given by 

(6.4) L {t, z) = g± {t, z)-^ Log± {t, z) , 
where g+ (t) , g- {t) solve the matrix Riemann problem 

exp t grad [Lq] (.z) = g+ {t, z) c/_ {t, 2)~\ 

(6.5) 

9+ it) e g+, g- it) e g^. 



Note that grad^n,m is a Laurent polynomial and hence is regular in the punc- 
tured Riemann sphere CPi\ ({0} U {oo}) = C\ {0}. Hence the factorization prob- 
lem (6.5) has the following geometric meaning. The projective line CPi is covered 
by two domains U+ = CPi\ {oo} and C/_ = CPi\ {0}. The function exptdh{L) is 
regular in [/+ n [/_ = C\ {0} and may be regarded as the transition function of a 
holomorphic vector bundle over CPi. Factorization problem (6.5) amounts to an 
analytic trivialization of this bundle. It is well known (see [12]) that not all vector 
bundles over CPi arc analytically trivial: each n-dimensional bundle breaks up 
into a sum of line bundles, and their degrees, di, d„ e Z form a full system of 
holomorphic invariants of the given bundle. In the language of transition functions 
this means that ex.ptdh{L) admits a factorization of the form 

exp tdh{L) = 3+ {z, t) d (z) g- [z, t)~^ , 
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where d{z) = diag (z''^, Thus formula (6.4) requires that all partial in- 

dices di, are zero. One can prove that this is true at least for t sufficiently 

small [13]. The exceptional values of t e C for which problem (6.5) does not admit 
a solution form a discrete set in C; at these points the solution L (t) has a pole: 
the trajectory of the Lax equation goes off to infinity. 

6.1. Riemann Problem and Spectral Curves 

Theorem 6.5 provides a link between the Hamiltonian scheme of Section 3 and 
the algebro-geometric methods of the finite-band integration theory (sec [14]). 
Namely, the formula (6.4) for the trajectories immediately implies that the Lax 
equations linearize on the Jacobian of the spectral curve associated with the Lax 
matrix. The proof is so short and simple that I would like to reproduce it here. 

Let g = Ql{n,C) and L{z) = J2-^i^^^-^i ^ 3^ ^ matrix- valued Laurent poly- 
nomial. Let us consider the algebraic curve Fq C C\ {0} x C defined by the 
characteristic equation 

(6.6) det{L{z) - A) = 0; 

we may regard z,X as meromorphic functions defined on Tq. Assume that the 
spectrum of L {z) is simple for generic z (this key technical assumption is satisfied 
in most applications). For each nonsingular point P cTq which is not a branching 
point of A there is a one-dimensional eigenspace E (P) C C" of L{z{P)) with 
eigenvalue A (P). This gives a holomorphic line bundle on Fq defined everywhere 
except for singular points and branching points. Let F be the nonsingular, compact 
model of Fq. One can show that the eigenvector bundle extends to a holomorphic 
line bundle E — > F on the whole smooth curve F.^ 

The spectral curve F with two distinguished meromorphic functions z and A 
and the line bundle E — > F constitute the set of spectral data for L{z). The 
evolution determined by a Lax equation of motion leaves the spectral curve (7.3) 
invariant, and the dynamics of the line bundle E is easy to describe. Let h be the 
Hamiltonian of our Lax equation; set = dh {L)\ since [L, M] = pointwise and 
the spectrum of L is simple, the eigenvectors of L [z) are also the eigenvectors of 
M{z), 

M{z {P))v = fi{P)v, P eT, V e E{P) c C", 

where fJ,{P) is a meromorphic function on F. Define the domains V± C F by 

T4 = |P e F;z(P)=^^ ^ oo|. Clearly, V+ U V- = F and ^ is regular in the 

intersection V+ HV- = Tq. Recall that a line bundle over a curve is specified by 
its transition function (with values in the multiplicative group C*) with respect 
to some covering; tensor product of bundles corresponds to the ordinary product 

^Indeed, the mapping P E (P) determines a meromorphic mapping V CPn-i (this is 
a corollary of the elementary analytic perturbation theory). By a standard theorem, any such 
mapping is actually holomorphic, and hence the eigenvector bundle extends to F. 
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of transition functions. The equivalence classes of line bundles form an abelian 
group PicF with respect to the; tensor product. Let be the line bundle on 
r determined by the transition function exp tji (P) with respect to the covering 
{y+ , V- } For alH e C the bundles Ft have degree zero and form a 1-parameter 
subgroup in the Picard group FIcq F of equivalence classes of holomorphic line 
bundles of degree zero on F which, by Abel's theorem, is canonically isomorphic 
to JacF, the Jacobian of F. (See, for instance, [15].) 

Theorem 6.6. The line bundle E regarded as a point o/PicF evolves linearly 

with time, E {t) = E ® F_t 

Proof. Since the Lax matrix evolves by similarity transformation, its eigenvectors 
evolve linearly. Let g± {t) be the solution of the Riemann problem (6.5). In view of 
(6.4), the moving eigenspace E^ (P) regarded as a subspace of F x C", is expressed 
as 

(6.7) EtiP)^g+{t,ziP))-'EiP) 
over V+, and 

(6.8) EtiP)=g_it,z{P))-'EiP) 

over In other words, g± [t, z) ^ define isomorphisms between E (P) and Ei (P) 
over V±. The transition function in V+nl^_ which matches these two isomorphisms 
is 5+ {t, z (P))~^ g- {t, z (P)) |_Et(p). It is easy to check that 

5+ {t,z)g- {t,zy^ = exptdh{LQ) 

implies 

g- {t, z) g+ {t, zy'^ = exptdh {L (t)) ; 

hence, 

(6.9) .9+(t,z(P))-i5-(t,^(P))U.(P) =e-*^(*'^(^»U.(P) =e-*'^(^), 

or Et (P) = £■ (P) (8> P_t, as claimed. □ 

The eigenvector ip {t, P) = (-01, V^„) G E{t,P) C C" is called the Baker- 
Akhiezer Junction of L{z,t). From (6.7, 6.8) it follows that the Baker-Akhiezer 
function ip {t, P) in the domains V± C F may be written in the form 

V'± {t,P)=g± {t,z{P))-^i^{P), 

so that 

(<,P) = e-*'^(^)V- {t.P). 
Since dtg± ■ g±^ = M± (see the proof of theorem 3.3), we have 



^ V± {t, P) = -M± {t, z (P)) V± {t, P) . 
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Using the machinery of Algebraic Geometry, it is possible to construct the Baker- 
Akliiczcr function explicitly, in terms of the Riemann theta functions and Abelian 
integrals. This, in turn, allows to obtain an explicit solution of the Riemann 
problem. Let us explain how the matrices g± {t, z) may be reconstructed from 

^± (^j P)- Suppose that z G C is not a ramification point of F (i.e., all eigenvalues 
of L{z) arc distinct); let Pi,...,P„ be the points of F which lie over z. Let us 
arrange the column vectors ij.)± {t. Pi) , ip± {t, P„) in a n x n matrix tp± {t, z). 
Put 

(6.10) 9±(i,^) =V'±(^,A)V'±(0,A)-\ 

Note that if wc change the ordering of branches Pi, P„, ^± is multiplied on the 
right by a permutation matrix and hence g± remains invariant. 

Proposition 6.7. (i) g± {t,z) satisfies the differential equation 

(6.11) ^^^^ = -M± {t, z) 9± {t, z) , M± {t, z) = (grad/z [L (i, z)])± . 
(a) g± are entire functions of z"^^. (Hi) g± solve the factorization problem 

9+ {t,z)g- {t,zy'^ = exptgrad/i[L(0,2;)] . 

The key assertion (ii) is an easy consequence of the differential equation (6.11); 
indeed, by (6.10) g± is the fimdamcntal solution of (6.11) normalized by g± (0, z) = 
Id; hence it is holomorphic in the domain where M± (t, z) is nonsingular, i.e., in 
CPi\{oo} and CPi\{0}, respectively. 

7. More Examples 



7.1. Twisted Loop Algebras 

The decomposition of the loop algebra we used so far is based on its standard 
grading (i.e., grading by the powers of the loop parameter z). This grading is by 
no means unique, and it is possible to use other gradings to produce new examples 
of Lax representations. Another possibility (which actually absorbs the former 
one) is to bring into play twisted loop algebras. Let a be an automorphism of q of 
order n. The twisted loop algebra L{q, a) is the subalgebra of Lq defined by 

(7.1) L{q, a) = {x& Lg; a {x (z)) = x (ez)} , 

where e = exp ^ is the root of unity. Equivalently, L{Q,a) C Lq is the stable 
subalgebra of the automorphism a : Lg -^Lg such that x"' {z) = a {x {ez)) for 
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all X G Lq. The isomorphism class of the twisted loop algebra depends on the 
properties of cr; one can show that when cr is an inner automorphism, the sta- 
ble subalgcbra L[g,a) is isomorphic to Lg as an abstract Lie algebra; however, 
the grading which is induced on L{Q,a) C Lq by the standard grading in Lq is 
different. The; classical theorem due to V.Kac asserts that all different gradings 
on Lq may be obtained in this way. (Among the integrable systems that may be 
very naturally constructed along these lines one may quote periodic Toda lattices 
which arc associated with the decomposition of loop algebras derived from their 
principal grading, see [8] for details.) 

The most interesting case for applications is when a is an outer automorphism 
of order 2 (an involution). In this case we may assume that g and Lq are real. 
Here is the key example: Q = Ql{n), a {X) = —X* ( t denotes transposition). The 
loop algebra Lq"^ consists of Laurent polynomials 

x(^) = ^Xfc^^ 

where X2p = —X2p, X2p+i — -'^2p+i- Antisymmetric matrices belong to the Lie 
algebra so{n) of the orthogonal group, which describes kinematics of the rigid 
body; on the other hand, symmetric matrices are reminiscent of the quadratic 
form associated with kinetic; energy ( inertia tensor). Thus, the twisted loop 
algebra seems to be a good candidate to set up the stage for applications to the 
mechanics of the rigid body. Let us equip Lq'^ with the inner product 

(7.2) {X, y ) = -Res ^=0 -tr X (2) y (z) 

z 

which sets Lq'^ into duality with itself (mind the difference with (6.1): the factor 

makes the coupling respect parity; the minus sign makes the inner product 
positive on so{n)). Let r be the classical r-matrix associated with the standard 
decomposition Lq"^ = Lq'^+Lq'^, as in (6.2). Obviously, 

[LqIY (LqI)^ = ©fc<o0 • z\ (LqI)* :^ (LqI)^ = (Bk>o9 ■ zK 

Coadjoint orbits of ifl^ are direct products of orbits of ifl^ lying in (Lg^)"*" 

and orbits of Lq'L lying in (Lg^)^. Lax matrices describing the motion of the 
rigid body and related mechanical systems belong to the simplest ad*-invariant 
subspace Lq"_^ ^ consisting of matrices 

L {z) = az^^ + I + bz,l e so{n) ,a = a\b = &*. 

Exercise 7.1. (i) All monomials bz e Lq'^ are 1-point orbits of Lq'L. (ii) Coadjoint 
representation of Lq'^ in the subspace {az~^ + l} C (Lq'^)* factors through its 
finite-dimensional quotient Lq'^/z'^Lq'^. 

(iii) The quotient algebra Lq^/z'^Lq^ is isomorphic to the semidirect product 
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Qq = so (n) K sym{n) of the orthogonal algebra so (n) and the space of symmetric 
matrices (with zero Lie bracket). 

The associated Lie group Gq is the semidirect product of SO (n) and the additive 
group of the linear space sym{n) of symmetric n x n-matrices. Its coadjoint orbits 
are easy to describe; here is a simple example: 

Exercise 7.2. Fix a unit vector e € M" and let a = e®e be the rank one orthogonal 
projection operator onto M-e C M". Let T*S"~^ be the cotangent bundle of the 
sphere 6""-^ =SO{n)-eC R" reaUzed as the subbundle of S"-^ x R", 

j,*gn-l ^ l^j.^^^ g gn-1 X ^n. ^ q_| 

There is natural map tt : T*S^~^ — > 0-\ onto the coadjoint orbit of Gq passing 
through the monomial z~^a, 

IT : {x,p) I — > z~^x <Si X + p Ax. 

{tt is actually a double covering.) 
The Lax matrix associated with this orbit has the form* 

L (z) = z~^x (g) X + p A .X + z 6, X e 5*0 (n) , p e M", (p, .x) = 0; 

the associated phase space describes the point moving on a sphere. The constant 
matrix 6 e sym (n) does not affect kinematics, but is quite useful to produce 
interesting Hamiltonians. The simplest Hamiltonian is 

= - iRes2=otr-L(z)^ = i - i (6x,x) ; 
4 Z 2. A 

it describes the so called Neumann problem (point moving on a sphere in a 
quadratic potential). 

Exercise 7.3. Describe the set of commuting integrals of motion for the Neumann 
problem which are quadratic in momenta. 

We may generalize this example in various ways; a useful remark is that we 
may avoid a too detailed description of the coadjoint orbits. Instead, we may 
produce a map tt onto such orbit, or a union of orbits, which is compatible with 
the Poisson structure but need not be a bijection (and so possibly introduces 
some extra variables). This idea is implemented in the following statement. Set 
K = SO (n) ,t = so{n) and let T*K ~ if x 6 be the cotangent bundle (equipped 
with its standard Poisson bracket). Fix a € sym (n) and consider the mapping 

7r:T*K^ ( gg)* : (fc, p) ^ k {p + a) kr' 

**We denote by p A a; the n x n-matrix with entries PiXj — PjXi\ in a similar way, the entries 
oi X are xiXj . 
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Proposition 7.4. n is a Poisson mapping (i.e., maps canonical Poisson brackets 
in T*K onto Lie-Poisson brackets in (Qq)* ); its image is a union of coadjoint 
orbits of the semidirect product Gq = K x sym (n) . 

[In this statement we identified with Mat (n) as a hnear space and also used 
the inner product {X, Y) = —trXY to identify the dual space Qq with Mat (n). ] 
The cotangent bundle T*K is naturally interpreted as the phase space of a rigid 
body in M"; we get a family of Lax matrices parametrized by points ofT*K : 

(7.3) L (z) = z-^kak-^ + kpk'^ + zb; 

Hamiltonians which may be derived from (7.3) include the so called Manakov case 
of the motion of a free top in R" (for n = 3 this is the classical Euler top), or more 
generally, the Manakov top in a quadratic potential. 

Remark 7.5. One may wonder, what is the relation of the low-dimensional 
Neumann system to the "big" phase space T*K (we have dimT*5"~^ = 2n — 2 
and dim T*ii' = n (n — 1)). The answer is that, for special choices of a S sym (n), 
the Hamiltonians associated with the Lax matrix (7.3) possess high symmetry 
(resulting from the redundancy introduced by tt). The Neumann system is the 
result of Hamiltonian reduction of the "big system" with respect to this symmetry 
group. 

One may notice that the use of the twisted loop algebra was indeed crucial: the 
built-in symmetry of the Lax matrix accounts both for the correct kinematics of the 
rigid body (antisymmetric matrices) and for the symmetry of the related quadratic 
forms (notably, of the kinetic energy). Further generalization is straightforward: 
we must scan the list of semisimple Lie algebras and their involutions and look for 
nice-looking opportunities. In this way we get the following list: 



Algebra 


Involution 


Related systems 


Ql{n) 


X ^ -X' 


Manakov top, Neumann system, etc. 


so{n, 1) 


X ^ -X* 


Lagrange top, spherical pendulum, etc. 


so {p, q),p>q>2 


X ^ -X* 


Kowalevski top and its generalizations 


so{n, n) 


X ^ -X* 


Interacting Manakov tops 


G2 C so(4, 3) 


X ^ -X* 


Exotic integrable top on 5*0(4) 



In all cases, Lax matrices belong to the subspace Lq'Li i, to get particular exam- 
ples (for instance, the Kowalevski top) one sometimes has to perform additional 
Hamiltonian reduction; wc refer the reader to [8] for details. 

Let us finally discuss the implications of the twisting automorphism for the 
geometry of the spectral curve and for the linearization theorem. 



Proposition 7.6. (i) Let us assume that the Lax matrix L{z) G Lq'^ and a 
is an inner automorphism o/fl, order = m. In that case the spectral curve T = 
{(-2, A) e C^; det (L (^) — A) = 0} admits an automorphism a : {z,X) 1— > {ez,X) 
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(here e = exp27r?/n is the root of unity); this automorphism lifts to PicF and the 
transition function (6.9) which determines the evolution of the eigenhundle of L is 
invariant under a. (ii) Suppose that a is an outer involution, a {X) = —X*. Then 
the spectral curve admits an automorphism a : A) i— > (— ^, —A); the transition 
function is anti-invariant under a, i.e., a : exp/:;U,(z, A) i-^ exp {—tii{z, X)). 

The check of both assertions is obvious: An inner automorphism preserves the 
eigenvalues of a matrix; by contrast, a (X) = —X* changes the sign of the eigen- 
values. The logarithm of the transition function n {z, A) is the eigonvahic of the 
gradient gradi/ [i] (z) (more precisely, one of its branches associated with the 
eigenvector of L{z) which corresponds to A). Since gradiJ [L] e Lq'^ , iJ,{z,X) is 
invariant when a is inner and changes sign when it's derived from transposition. 

7.2. Rational Lelx Matrices 

Let us now describe how to deal with Lax matrices which are rational functions 
of the spectral parameter. As we mentioned, it is possible to trace down the 
corresponding coadjoint orbits inside the dual space of an appropriate completion 
of the standard loop algebra, but it is more practical to choose our basic Lie algebra 
in a different way. 

Let D = {zi, zn} C CPi be a finite set; we assume that oo G D. For Zj G D 
let Aj be the local parameter on CPi at 2: = Zj, i.e., Xj = z — zj if zj ^ 00 
and Aoo = z^^ for zj = 00. We define the local algebra 25^. as the algebra of 
formal Laurent series in local parameter with coefficients in a little Lie algebra 
Q, (Szj = g {{Xj)). (We may assume that g is the matrix algebra with the standard 
inner product.) If Zj ^ 00, let 25+ be the algebra of formal Taylor series in local 
parameter; for zj = 00 we set 25+ = Aoofl [[Aoo]] (in other words, 0+ consists of 
formal Taylor series without constant term). Put 

(direct sum of Lie algebras). Let q{D) be the algebra of rational functions on CPi 
with coefficients in g which are regular outside D; it is naturally embedded into 
©D (the embedding assigns to each X e g{D) the collection of its Laurent series 
at each point of D). 

Proposition 7.7. (i) (Sd = 9iD)+(3~jj (direct sum of linear spaces), 
(ii) The C-bilinear inner product on *8d 

(7.4) {X,Y) = ReS;,jtT XjYjdXj 

zjeD 

is invariant and nondegenerate. (Hi) g{D) and 25^ are isotropic subspaces with 
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respect to (7.4); moreover, 0(1))— (25^) . (iv) Coadjoint orbits of in q{D) are 
finite- dimensional. 

Sketch of a proof. An clement X = (Xj)^ is a finite collection of Laurent 
series; stripping each of them of its positive part we get a set of principal parts 
at Zj e D; let X° be the unique rational function with these principal parts; by 
construction, X — X^ £ (3~jj (mind the special role of cxd which fixes the normal- 
ization condition!). In brief, we can say that the decomposition (8d — q{D)+<3^ 
is equivalent to the Mittag-Leffler theorem for rational functions. Isotropy of g(D) 
and ©J means that the inner product restricted to these subspaces is identically 
zero; this condition assures that the associated classical r-matrix is skew symmet- 
ric. For (8~jj this assertion is immediate, since (for z ^ oo) the product of two 
Taylor series has zero residue; for z = oo the residue disappc^ars because of the 
normalization condition. The isotropy of g{D) is a reformulation of the classical 
theorem: the sum of residues of a rational function is zero. 

Since ®^ is a direct sum of local algebras, its coadjoint orbits are direct products 
of the coadjoint orbits of each local factor; it is easy to see that the coadjoint orbits 
of the local algebra 0^. are modelled on rational functions with a single pole at 
z = Zj] moreover, the subspace of rational functions with prescribed order of 
singularity at this point is stable under the coadjoint action of 25+ . Clearly, this 
subspace is finite-dimensional, which proves (iv). 

Exercise 7.8. Describe coadjoint orbits in the subspace of functions admitting 
only simple poles. 

Our main theorem immediately applies in this setting and provides an ample set 
of integrals of motion in involution for Lax equations with rational Lax matrix. 

Remark 7.9. (1) Since we are interested only in Lax operators which are global 
rational functions on the Riemann sphere, we consider only coadjoint orbits of 
25^ C (25d)^; this is legitimate, since (25 zj)^ splits into direct sum of two comple- 
mentary subalgebras, 

(25z,),~6+ ©0(D), 

and hence its orbits are direct products of orbits lying in g{D) and in &~^; in other 
words, we take orbits which project into zero in 25 J. 

(2) The global algebra 25^ is decomposed into direct sum of local factors, 
ffizje£>25+; coadjoint orbits of each local algebra are the same that we encoun- 
tered for the ordinary loop algebra. What makes things different, is the way these 
orbits are embedded into the bigger algebra; this embedding affects the choice of 
the invariant Hamiltonians as well as the formulation of the factorization problem. 

The use of formal series is well adapted for the study of coadjoint orbits in g{D); 
in order to be able to define Lie groups associated with our Lie algebras, we must 
change the topology by replacing formal series with convergent ones. Let 
be the group of germs of functions with values in G which are regular in some 
punctured disc around Zj G CPi (with topology of uniform absolute convergence), 



32 



Semenov-Tian-Shansky 



{Qtj)^ ^ QYj its subgroup consisting of functions regular in the entire small 
disc, and G{D) the group consisting of holomorphic mappings CPi\Z) G. The 
infinitesimal decomposition described in proposition 7.7 (i) corresponds to the 
following multiplicative problem: 

Given a set of local meromorphic functions g\, ■■■,gN, gj € G^, find a global 
meromorphic function go which is regular in the punctured sphere CPi\D such 
that gogj^ is regular in some small disc around Zj. 

This is the standard multiplicative Cousin problem; its geometrical meaning is 
the same as for the matrix Riemann problem discussed above: it corresponds to 
the trivialization of a vector bundle over CPi (defined with respect to a different 
covering of the sphere). 

Exercise 7.10. Reformulate the global factorization theorem in this setting. 

One is of course tempted to generalize the above construction replacing CPi 
with an arbitrary Riemann surface. There is an obvious obstruction which comes 
from the Mittag-Leffler theorem for curves: a global meromorphic function on a 
curve r with prescribed principal parts exists if and only if these principal parts 
satisfy a set of linear constraints; roughly, the sum of residues 

Res^^ trXjo; 

must be zero for all holomorphic differentials w e -ff^(r) ® Q- The trouble is that 
the constrained data do not form a Lie subalgebra inside the global algebra ©d, 
and hence one cannot find a complement of which is a Lie subalgebra. When 
r is elliptic, this obstruction may be overcome by imposing additional automorphy 
conditions, see e.g. [8]. 



8. Zero Curvature Equations 

In applications to integrable PDE's, Lax matrices are replaced by first order ma- 
trix differential operators. The systematic treatment of these applications is based 
on the use of double loop algebras^ or, more precisely, of their central extensions. 
Let us start with discussion of the central extension of the ordinary loop algebra. 

8.1. Central Extensions of Loop Algebras 

Set Q = fll(n) and let = C°° {S^;q) be the Lie algebra of smooth functions 
on the circle with values in q and with pointwise commutator. (Mind that in the 
present setting we choose topology in our loop algebra in a different way! This is 
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because we are willing to treat functions of x G as dynamical variables for our 
future evolution equations.) We equip & with the invariant inner product 



(8.1) {X,Y)=[ 

Jo 



2f! 

irXYdx: 



accordingly, we get an embedding (S C 0* which defines the smooth dual of (S. 
Put 

/■27r .y 

(8.2) w{X,Y)= / tvX-—dx; 

Jo dx 

Exercise 8.1. w is a skew symmetric bilinear form which satisfies the cocycle 
condition 

(8.3) uj {[X, Y],Z)+uj {[Y, Z],X)+uj {[Z, X],Y) = 0. 
Put = 65 e C and define the Lie bracket in S by 

[{X,c),{Y,c')] = {[X,Y],u{X,Y)), x,Yee, cc'gC. 

(The Jacobi identity is equivalent to (8.3).) Notice that c = {(0, c) ; c G C} C © 
is central in 6 and © may be identified with the quotient algebra, 6 = ©/c. The 
(smooth) dual of & may be identified with (S ® M. 

Proposition 8.2. The coadjoint representation of ® is given by 

(8.4) ad*X ■{L,e) = {[X,L] + ed^X,0). 

(Note that this representation is trivial on the center c C © and therefore it may 
be regarded as a representation of&.) 

Let G = C°° [S^; G) bo the Lie group associated with the Lie algebra ©. The 
coadjoint representation of © in ©* may easily be integrated to a representation 
of G. 

Proposition 8.3. The coadjoint representation of G in ©* is given by 

(8.5) Ad*g-{L,e)={gLg-^ + ed,g-g-\e). 

Proposition 8.3 admits a remarkable geometric interpretation. Consider the 
auxiliary linear differential equation 

(8.6, e^ = L^ 

(we regard it as a differential equation on the line with periodic potential L). 
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Exercise 8.4. Coadjoint representation (8.5) leaves invariant the hyperplanes 
e = const in &*; on each hyperplane e = const ^ it is equivalent to the gauge 
transformation of the potential L in the linear equation (8.6) induced by the natu- 
ral action of G = C°° (S^; G) on its solutions by left multiphcation, g : tp g ■ ip. 

Let tpo € G°° (M; G) be the fundamental solution of (8.6) normalized by tjjo (0) = 
id; Tl = tpo (27r) G G is called the monodromy matrix of L. 

Theorem 8.5. (Floquet) Two potentials L,L' e C°° (S^;q) lie on the same 
coadjoint orbit in &* (with fixed e ^ 0) if and only Tj, and T^i are conjugate in 
G. 

Sketch of a proof. Without restricting the generality we may assume that Tj, = 
T^'. Let V'-LjV'l' be the fundamental solutions normalized by tpj^ (0) = tpL' (0) = 
id; put g{x) =tpL {x) V'l' Clearly, g is 27r-periodic on the line and g ■ ipL' = 

ipL, hence Ad*g ■ L' = L. 



Corollary 8.6. All coadjoint orbits lying in the hyperplanes e = const =^ have 
finite codimension; the ring of Gasimir functions is generated by spectral invariants 
of the monodromy. 

The Hamiltonian mechanics in 25* may be defined with the help of the elementary 
calculus of variations. Let [L] be a smooth functional of the potential L, = 
gradt^ [L] e (S its Prechet derivative defined by 



d f^'^ 
—tp[L + sri]= tr (x) rj (x) dx. 
ds Jo 



The Lie-Poisson bracket of two functionals ipi , ip2 is given by 

{Vi, ¥'2} [L] = tr (L (a;) [X^^ (x) , X^^ (x)] + eX^^d^X^^) dx. 
Jo 

Proposition 8.7. The Hamiltonian equation of motion on &* with Hamiltonian 
(f is equivalent to the following differential equation for the potential L : 



dt dx 

Equation (8.7) has a nice geometrical meaning. Let us consider the g- valued 
differential form 

Ldx + X^dt; 

it may be regarded as a connection form of a connection on (with values in 
g); equation (8.7) then means that this connection has zero curvature (hence the 
term "zero curvature equation" .) We would like to use the central extension 6 as 
a building block to construct integrable equation; there are already two reassuring 
points: 
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1. The description of coadjoint orbits in (S* automatically leads to the auxiliary 
linear problem (8.6). 

2. Equations of motion associated with (S are zero curvature equations, as de- 



However, there are also two major drawbacks: 

1. There is only a finite number of independent Casimirs (one can take e.g. the 
coefficients of the characteristic polynomial det {Tl — A)). 

2. The Casimirs are highly nonlocal functionals of the potential. 

By contrast, in order to get integrable PDE's we need an infinite number of con- 
servation laws; these conservation laws are usually expressed as integrals of local 
densities which arc polynomial in the matrix coefficients of the potential L and its 
derivatives in x. The way to resolve these difficulties is suggested by the auxiliary 
linear equation (8.6): in order to characterize the potential, we need to know the 
monodromy /or all energies: in other words, we must introduce into (8.6) spectral 
parameter. Algebraically, this means that we have to modify the choice of our 
basic Lie algebra. 

8.2. Double Loop Algebras 

Let us put © = C°° (S'^fl) as before; let g =(S(g)C{{z)) be the algebra of formal 
Laurent scries with values in ©. (In other words we can say that g is the double 
loop algebra of q; the different roles of the variables x, z are imposed by our choice 
of its central extension.) We equip g with the inner product 



(8.8X^,1^) =Res^=o / tr X{x, z)Y {x,z)dx = — / ti X {x, z)Y {x,z)dxdz. 



Let g be the central extension of g defined by this cocycle.® As before, we may 
identify the dual of g with g ® C. Formula (8.4) for the coadjoint representation 

remains valid. We conclude that the coadjoint representation for g coincides with 
(infinitesimal) gauge transformation associated with the auxiliary linear problem 

®In the case of the loop algebra the definition of the 2-cocycle uj is essentially unique; for the 
double loop algebra there is a possibility to introduce into (8.9) a weight factor (j>{z) which does 
not depend on x; this weight factor will modify the auxiliary linear problem. This freedom is 
useful in applications; our choice in (8.9) is the simplest one possible. 



sired. 




The 2-cocycle w on g is defined by 



(8.9) 
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(8.6), where this time L G g, i.e., it is a formal series in z with coefficients in 
C5. In other words, z plays the role of spectral parameter in the auxiliary linear 
problem, as desired. There are some troubles with the definition of the associated 
Lie group, but let us ignore them for the moment. Notice that if i is a polynomial 
in z, z~^, the monodromy matrix Tl is a well-defined analytic function of z (with 
values in G = GL (n)) which is holomorphic in C\ {0}. 

Our choice of the basic Lie algebra makes it easy to define the other key element 
of our scheme, the classical r-matrix. Set g+ = 25 ® C [[z]] , g_ = (S ® z~^C [z~^] ■ 
Clearly, g = g+-i-g_ as a linear space. Both subalgebras g+ and g_ are isotropic 
with respect to the inner product (8.8) which sets them into duality. As before, 
wc put 



where P+,P_ are the projection operators onto g+ and g_, respectively, and 
define the r-bracket on g by [-'f, = 5 ([''-'^i^] + i-'^)'''^])- In this way we get 
the algebra g^, but this is still not quite what wc need to apply our main theorem, 
since our basic algebra is g, not g! As it happens, the theorem survives central 
extension. 

Lemma 8.8. Let u) be a 2-cocycle on g and r € Endg a linear operator which 
satisfies the modified Yang-Baxter identity. Set 



Then ujr is a 2-cocycle on gr. 

Exercise 8.9. Prove lemma 8.8 (the proof is abstract and uses only manipulation 
with the Jacobi and the Yang-Baxter identities). 

Let gr be the central extension of gr associated with ujr- It is easy to see that 
(g, gj^) is a double Lie algebra and wc may apply our main idea: Casimirs of g 
regarded as Hamiltonians with respect to the Lie-Poisson bracket of g^ give rise 
to generalized Lax equations. Actually, there is one more simplification, which is 
due to our choice of w (see (8.9)): 

Exercise 8.10. Suppose that r G Endg is skew symmetric with respect to the 
inner product on g; then = 0. 

The r-matrix (8.10) clearly satisfies this condition; hence the algebra g^ = gr-©C 
splits and the Lic-Poisson brackets for gP and g^ coincide. Explicitly, this means 
that the Poisson bracket of two smooth functional ipi , <^2 defined on g* ~ g is 
given by 



(8.11) {(pi,(f2}rm = — / tr ([grad(pi,grad(y9i]^(x,z) • L (x, z)) dxdz. 



(8.10) 



r = P+ 



Ur {X, Y) = ^{oj {rX, Y)+u {X, rV)) . 




where gradi^j [L] {x, z) G g, i = 1, 2, is the Prechet derivative. 
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RerriEirk 8.11. The skew symmetry of r makes the above discussion of the cocycle 
ojr void; however, as aheady noticed (see footnote 9), we may modify the cocycle 
a; by a weight factor (/)(^;), and in that case our Poisson bracket will contain 
derivatives of the gradients. 

The antisymmetry of r makes the description of coadjoint orbits very simple. In 
the absence of cocycle we must deal with the orbits of gr; note that the r-matrix 
(8.10) is acting only on the variable z in the double loop algebra, and hence the 
other variable x becomes a parameter. Let us consider the 'little' loop algebra 
ifl = (S) C {{z)) and the associated algebra Lq^ which we have already discussed 
in Section 6. The 'big' algebra consists of smooth periodic functions on the line 
with values in Lq,. and with pointwise commutator. 

Proposition 8.12. Fix an orbit O C (Lgr)*/ then O = Map{S'^,0) is an orbit 

Ofgr- 

More generally, we may vary orbits O lying over different points of (i.e., 
consider smooth families of orbits in (iflr)* parametrized by S^). 

Example 8.13. Let g = si (2); then the matrices s G fl, 



form a coadjoint orbit <Si C sl{2). Check that Oh = s,s G Si} c ~ 

Lg* is a coadjoint orbit of Lg+ C -Lflr- The corresponding orbit Oh C g is 
parametrized by a triple of 27r-periodic functions Sj,j = 1,2,3, satisfying the 
constraint (8.12). The associated linear differential operator is 



One can show that the simplest local Hamiltonian associated with (8.13) is 



the corresponding nonlinear equation describes the Heisenberg ferromagnet: 



(8.12) 





— [^i ^xx] • 



Example 8.14. Set 




-"^'^ Poisson bracket of functionals which docs not contain derivatives dx of the gradients is 
sometimes called ultralocal; in more complicated cases, Poisson brackets may contain derivatives 
(non-ultralocal case) or even be non-local, i.e., contain integral operators. 
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The matrices 

^ V 

form a coadjoint orbit of the Lie algebra Lq- C LQr- The corresponding orbit 
Os C g is parametrized by a pair of functions u,v € C°°(S'^,C); the associated 
linear operator is 

dx ^' 



this is (essentially) the Lax operator for the so called split nonlinear Schroedinger 

idtu = -Uxx + u^v, - idtv = -Vxx + v^u. 



equation^ ^ 



In a more general way, for q = Qi{n) we can easily construct coadjoint orbits 
associated with linear differential operators of the form 

(8.14) j.-u{x,z), 
where the potential is a Laurent polynomial in z. 



8.3. Monodromy Map 



Examples 8.13, 8.14 show that our approach is indeed reasonable: starting with 
double loop algebras, we arrived at the natural auxiliary linear problems with spec- 
tral parameter of the form which is familiar in spectral theory. The next question 
is to construct an appropriate class of Hamiltonians. Formally, the Hamiltonians 
are spectral invariants of the auxiliary linear operator, i.e., the spectral invariants 
of its monodromy matrix. The monodromy matrix Mu of (8.14) is a holomorphic 
function in CPi\ ({0} U {oo}); any functional of the monodromy which is invariant 
under conjugation gives rise to a zero curvature equation on the coadjoint orbits 
of gr. The mapping M : [/ ^ Mjj is the direct spectral transform for (8.14). Wc 
may regard M as a mapping from g into the group of matrix-valued functions 
which are regular in the punctured Riemann sphere. When the Poisson structure 
is ultralocal (i.e., the r-matrix is skew with respect to the inner product in g), 
the spectral transform M has an important property: the target space carries a 
natural Poisson bracket and M is a Poisson mapping (i.e., it preserves Poisson 
brackets). Wc shall briefly outline this computation, since it plays an important 
role in the theory and brings into play an important class of Poisson brackets (the 
so called Sklyanin brackets). To simplify the notation, let us start with ordinary 
loop al,t!,'ol)ras. 

^-'^For simplicity, in both examples we deal with complex Lie algebras; in order to pass to a real 
form we must choose an anti-involution of our basic loop algebra; this is of course a necessary 
step in order to make the auxiliary linear operators genuine skew selfadjoint operators. 
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Let g = Ql{n), and r S End 5 a linear operator which satisfies the modified Yang- 
Baxter identity (4.4) and is skew symmetric with respect to the inner product on 
0. Let (S = C°°{S^,q); we equip 6 with the inner product (8.1); the Poisson 
bracket of functional on (S* ~ (S is given by 

/•27r 

(8.15) {</9i,</92}r [L] = / tr ([grad</9i [L] (a;) ,gradv?2 [L] {x)]^L{x))dx. 

Jo 

For i e (S let tpL be the fundamental solution of (8.6) normalized by i/'l(0) = id 
and M]^ G GL{n) the monodromy matrix. Fix a smooth function ip € C°° {GL{n)) 
and put [L] = (f {Ml). Let V<^, VV € be the left and right gradients of on 
G = GL (n) defined by 

4- p{e'^x) = tr (X • Vip {x)), ^ ip{xe'^) = tr (X • VV (x)) ,X GQ. 
ds s=o as «=o 

Lemma 8.15. T/ie Frechet derivative of the functional h^p is given by 

(8.16) grad/i^ [L] {x) = (a;) VV (Ml) ^ (x)"' . 

CoroUciry 8.16. The Frechet derivative satisfies the differential equation 

(8.17) f^i^,^] 
and the boundary conditions 

(8.18) X (0) = W (Mi) , X (27r) = (Mj,) . 

Sketch of a proof. Taking variation of the both sides of (8.6), we get 

(8.19) dML=L5^JL + {SL)^JL- 

Set 5'4)L = i^hY, where Y is an unknown function, Y e C°° (M, g) , F (0) = 0; 
(8.19) yields dxY = il)~^6Ltp, whence 

Y{x)= r i^l'{y)SL{y)^L{y)dy. 
Jo 

Since Ml = tpL (27r), we get 

M£'SMl= [ iJl' -SL-i^Ldy. 
Jo 

Now, 

S>fi {Ml) = tr (W {Ml) ■ M£^6Ml) = [ tr V'L(t/)VV {Ml) i^l\y) ■ SL{y) dy, 

Jo 
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which impUes (8.16). Taking derivatives of the both sides of (8.16) yields (8.17). 

Proposition 8.17. The Poisson bracket of two functional h^^,h^2 on C5* is 
given by 

(8.20) {V,/i^J[L] = /i{^„^,}[L], 
where the Poisson bracket of , <^2 S C°° (G) is defined by 

(8.21) Wu^2}g = ^tr (r (Vtpi) - r ( Wi) VV2) • 

Corollary 8.18. Let us equip the group G = GL (n) with the Poisson bracket 
(8.21). Then the monodromy map M : (S* — > G : L Ml preserves Poisson 

brackets. 

Proof. Set Xi = grad/i,^^, i = 1,2. We have 

/.2-n- 

{V,/i^J[L]=/ tY {[Xi,X2lL)dx 
2 

= ^J^ tr {[rXi,X2] + [Xi,rX2])Ldx 
= I r tr ([L, X2] rXr + [L, X^] rX2) dx 
-tr {rXi ■ X2) dx, 



2/ 



d 

dx ' 



where we used the definition of the r-bracket, the cyclic invariance of trace, the 

differential equation (8.17) satisfied by Xi and, finally, the skew symmetry of r. 
Evaluating the last integral and taking into account the boundary conditions (8.18) 
for Xi, we get (8.20). □ 

Formula (8.21) defines a remarkable Poisson bracket {the Sklyanin bracket) on 
the group manifold itself. The Jacobi identity for this bracket is not obvious 
(though it follows from our computation). Its properties will be discussed in some 
detail in Section 9. 



Exercise 8.19. Show that central functions on G (i.e., functions which satisfy 
if (xy) = (p {yx) for all x,y G G) commute with respect to the Sklyanin bracket. 

The above discussion applies to loop algebra 05 — C°°{S^,g); the generaliza- 
tion to the double loop algebra is straightforward: in our computation, we must 
replace the finite dimensional algebra g with its loop algebra Lq; accordingly, 
smooth fimctions ipi,ip2 G (G) are replaced by smooth functionals on the 
corresponding loop group, their left and right gradients are replaced by the left 
and right variational derivatives, etc. Spectral invariants of the auxiliary linear 
problem (8.14) correspond to central functionals on the loop group. An example 
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of such a functional is given by evaluation functionals Hn w [U] = tvMJ}{w), w € 
CPi\({0}u{oo}). 

Exercise 8.20. Compute the variational derivative of ifn,^ with respect to U. 

The fundamental drawback of these functionals is, however, their nonlocality. 
The description of local functionals is outlined in the next paragraph. 

8.4. Formal Diagonalization and Local Conservation Laws 

In contrast to evaluation functionals, local conservation laws arc related to the 
asymptotic expansions of the monodromy matrix at its essential singularities, i.e., 
for 2 = 0, 00. This implies some additional difficulties: 

1. These functionals arc not defined everywhere on the double loop algebra. 

2. The associated formal scries arc in most cases divergent. 

Let us assume that the potential U{x,z) in the auxiliary linear problem (8.14) 
is a Laurent polynomial, 

M 

U = Y,Ukz\Uk&C'^ {S\Qi{n)); 

-N 

let Jo = U-M, Joo = Um be its lowest and highest coefhcients.. 

Definition 8.21. U is called regular if 

(i) the matrices Jq {x) , Jqo (x) are semisimple, 

(ii) the centralizers of Jo (a;), Joo{x) in g = Ql{n) are conjugate for all x G S^. 

We have seen that the set of Laurent polynomials of fixed degree is a Poisson 
subspace for the r-bracket. It is easy to check that the regularity condition holds 
for entire coadjoint orbits in this subspace, and hence it is a characteristic of our 
phase space. Our next theorem allows to construct for regular Lax operators two 
series of local Hamiltonians which are associated with the poles of U on CPi . For 
concreteness, we shall describe the construction of the series associated with the 
pole at infinity. Performing a suitable gauge transformation we may assume that 
the leading coefficient at infinity Joo {x) satisfies a stronger condition: 

(W) the centralizer of Joo in q = Ql{n) is a fixed subalgebra C Q which 
does not depend on x. 

(Note that since, by construction, local Hamiltonians are gauge invariant, this 
stronger condition does not restrict generality.) Let qj^ C 0''°° be the commutant 
of 0"'- , 

QJoo ={Xg&; [X, Y] = for all Y e q^-} . 
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Theorem 8.22. (On normal form at infinity) There exists a formal gauge trans- 
formation 

oo 

= + 5^ ^mz-"", e {S\Mat (n)) , 

rn—l 

which transforms the differential operator dx — U into normal form, 
(8.22) ($°°)"^ o ( -^^ - [/ ) o$°° = 



dx J dx 
where 

oo 
m=-M 

matrix coefficients of^^, are expressed as polynomials of the coefficients of 
U and its derivatives in x. 

Sketch of a proof. The intertwining relation (8.22) is equivalent to the differential 
equation 

(8.23) (^-^ -U^^°° = 

which may be solved recurrently in powers of the local parameter z^^. The first 
nontrivial coefBcients $i,£>_m+i satisfy 

(8.24) Joo$r - '^TJoo = D°^M+i - Um-1. 

This equation for $i admits a solution if and only if the r.h.s. is in the image of 
ad J oo S Endg. Assumption (i) above implies that 

= Imad Joo+kerad Joo, 

moreover, by (ii') ImadJoo and keiadJoo = 3"^°° do not depend on x. Hence, 
e g^°° is uniquely determined from the solvability condition of (8.24) and 

$r = (ad Joo)"' (^-M+i - Um^i) ■ 

If the coefficients Z)^^,j_,_-^, D'^j^_^^ are already determined, we get 

for ^m+i the relation of the form 

(8.25) ad Joe • $^+1 = -Frr, {U, $r , ^-M+l, D^M+m) , 

where F^^ depends on U and on the already determined coefBcients and their 
derivatives. By the same argument, (8.25) allows to determine D^M+m+i^^m+i- 

^^In various applications Jo, Joo arc regular matrices with distinct eigenvalues. In that case 
g'^o = Bjq and 0-'°= = Qj^ arc abclian subalgcbras; hence theorem 8.22 means that the potential 
U may be transformed to diagonal form by a formal gauge transformation. When the eigenvalues 
of Jo , Joo have multiplicities, the potential may be transformed only to block diagonal form. 
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Remark 8.23. The coefBcients ^^,D^ are determined from (8.25) not com- 
pletely canonically, since we must fix in some way the operator {adJoo)~^ ■ One 
can show that this freedom corresponds to the possibility to perform gauge trans- 
formations 

(8.26) - exp (-(/.) o I - ) o exp (j), 

dx \dx I 



Id+Y^ (t>mz eg-'- 

m=l 



The formal series is sometimes called the formal Baker function at infinity 
of the operator L = — U. For a G Qj^ (g) C [z, z~^] and put 

(8.27) [U] = Res ^=o / tr a {z) D°° {x, z) dx. 

Jo 



Theorem 8.24. 

(i) Functionals do not depend on the freedom in the definition of the normal 
form. 

(ii) All functionals are in involution with respect to the Poisson bracket 
(8.15) ong;. 

(Hi) Hamiltonian equation of motion defined by on g* have the form of zero 
curvature equations. 



Lemma 8.25. Gauge transformations (8.26) leave the density ti a (z) D°°{x,z) 
invariant up to a total derivative. 

Sketch of a proof. Gauge transformations (8.26) map D°° into e~'^D°°e'^ — 
e~'^dxe'^. By a standard formula, 

g-ad<f> _ / \ \ \ 

dx (exp0) = — dxcj) = [id - -ad(j) + — (a#)^ + ... • dxcj). 



-adcj) y 2^ 3! 

Hence 

/g-adct>_j^ \ 

tr a exp (-0) dx (exp (j)) = tr a • dxf/)] = 

\ —aa(p J 

/ e~"'^'^ — Id \ 
tr ( — aj ■ 8x4) = tr a • dx(t) = dx (tr a^) ; 

where we also used the invariance of trace, the condition a G gj^ which assures 
that it commutes with 4> and, finally, the condition dxOi = 0. 
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Lemma 8.26. The Prechet derivative of is given by 

(8.28) gradil~ = ($°°)"^ , 

where is the formal Baker function. 
Sketch of a proof. Taking variations of both sides of (8.23), wo get 

= ($°°)-i Jf/$°° + £)°°,($°°)"^^$°° -a^ (^($°°)"^^$°°) . 

Hence 
SH^ = 

Res^=oy |tr$~a($~)"^5;7 + tr (5^a - a]) ($°°)"^ dx, 

where we used the invariance of trace and integrated by parts; the contribution of 
the second term vanishes, since dxCt = [D°°, a] = 0. □ 

Corollary 8.27. The Frechet derivative X = gvadH^ satisfies the differential 

equation 

(8.29) d^X = [U,X]. 
Indeed, (8.28) and (8.23) imply 

S^X = (C/$°° - $°°£)°°) a ($°°)"^ - ($°°)"^ (C/$°° - $°°£)°°) ($°°)-^ 

= [U,X] - [£'°°,Q!] ($°°)-^ = [C/,X] . 
Note that geometrically (8.29) is equivalent to 

(8.30) grad H^[U]-U = 0, 

where adX is the coadjoint representation of the Lie algebra g (the central exten- 
sion of g) ; this is precisely the property which characterizes the Casimirs of a Lie 
algebra (cf. proposition 2.6). In the present setting is not a true Casimir 
function: it is defined only for regular elements U € g with fixed highest coeffi- 
cient. However, a short inspection of the proof of theorem 3.1 shows that it uses 
only (8.30); the last assertion of theorem 8.24 now follows. 

In a similar way, we may define the second series of local Hamiltonians which 
is associated with the pole at 2; = 0; one can show that the Hamiltonians from 
these two families mutually commute (this does not follow immediately from the 
arguments above, but may be proved in a similar way). In a more general way, 
if the potential U is & rational function on CPi(cf. Section 7.2), we may asso- 
ciate a series of local Hamiltonians to each of its poles; the corresponding Frechet 
derivatives are formal Laurent series in local parameter at the pole. 
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8.5. Higher order differential operators 



In applications, it is quite common to deal with Lax representations which con- 
tain higher order differential operators; the most famous example is the KdV 
equation associated with the Schroedinger operator on the line 

In order to put these operators into our framework, one needs extra work. We 
shall outline the procedure without going into details. First of all, an n-th order 
differential equation 



(8.31) 



+ ... + uoip + zip = 



may be written as a first order matrix equation, 

d 



(8.32) 
where 

(8.33) 



dx 



( 



+ = 0, 



\ 







1 





/ V 






i,' 

















V Mo + -2 

However, the companion matrix in (8.33) contains "too much zeros" and cannot be 
directly associated with a coadjoint orbit of the loop algebra. Let us observe first of 
all that choosing the column vector in (8.33) in this particular form is not quite 
canonical: we may add to fc = l,2,..., n — 1, an arbitrary linear combination 
(possibly, with variable coefficients) of i/i, ?/;', ^/j^^"-'^); this freedom amounts to a 
gauge transformation (p{x) ^ n {x)-ip {x), where n is a lower triangular (unipotent) 
matrix. The potential L in (8.33) becomes an arbitrary matrix of the form 



/ 



(8.34) 



L = 



1 



\ 







1 

* / 



\ * • • • * 

The companion matrix in (8.33) is the result of gauge fixing; indeed, we have 



Proposition 8.28. For each potential L of the form (8.34) there exists a unique 
lower triangular gauge transformation which transforms it into the canonical form 
(8.33). 
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Conclusion: The space of n-th order differential operators is the quotient space 
of the set of all potentials L of the form (8.34) modulo the ga.uge action of the 
lower triangular group. This quotient space is modelled on the set of companion 
matrices. With a little skill in symplectic geometry one may describe this quotient 
space in terms of Hamiltonian reduction (the key point is to observe that potentials 
of the form (8.34) form a level surface for the moment map associated with our 
gauge action). 

There is one more difficulty: the term of highest degree in z in the potential L (z) 
is a nilpotent matrix, and so the expansion procedure which yields local integrals 
of motion does not work. To tackle with this trouble let us recall that the loop 
parameter z is in fact associated with a special grading (the standard grading) of 
the loop algebra; in this grading, the constant matrix in (8.32) is 



(8.35) 



/O 1 












\ 





1 

) 



If wc use the so called principal grading of the loop algebra instead, (8.35) is 
replaced, after rescaling, with 



/ 





V 1 







\ 





1 

/ 



which is already a semisimple matrix with different eigenvalues; now local conser- 
vation laws may constructed by expansion in in the usual way. 



8.6. Dressing and Solutions of Zero Curvature Equations 



Formula (8.28) for the Frechet derivative of a local Hamiltonian makes it clear 
that the direct use of the global factorization theorem (theorem 3.3) to solve zero 
curvature equations is impossible. Indeed, gradiJ^ [U] is a formal series in local 
parameter z~^] it is therefore impossible even to define the 1-parameter subgroup 
exptgradiJ^ [C/]; this reflects real analytic difficulties which exist in the study 
of the initial value problem with arbitrary initial data for integrable PDE's. Our 
way around this difficulty is to introduce the following definition. 
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Definition 8.29. A difiFerential operator dx — U is called strongly regular at zero 
(at infinity) if U satisfies the conditions imposed in definition 8.21 and, moreover, 
its formal Baker function at zero (at infinity) is convergent. 

One may of course strongly doubt the merits of this definition. The point is, how- 
ever, that strongly regular operators form a homogeneous space: there is a natural 
action of the loop group (called dressing transformations) on the set of potential 
which preserves strong regularity. Moreover, commuting Hamiltonian flows gen- 
erated by local Hamiltonians are naturally included into the dressing group as a 
maximal commutative subgroup. Examples of strongly regular operators include 
solitonic and finite-band solutions; on the other hand, it is very difficult to give 
their complete characteristics in local terms (i.e., to tell which initial data in the 
space C°°(>S'^, (S) correspond to strongly regular operators). The group action in 
question was discovered by Zakharov and Shabat [16] (though at the time they did 
not notice the composition law), and later by Sato and his school [17]. A subtle 
question in the theory of dressing transformations is the treatment of boundary 
conditions. In our discussion above we used periodic boundary conditions; the 
main motivation was the Floquet theorem which gives an accurate description of 
coadjoint orbits and Casimir functions for the central extension of the loop alge- 
bra. Dressing transformations do not preserve periodicity.^^ One more interesting 
point is the relation between dressing and the Poisson brackets: since the dressing 
action is defined in terms of the Riemann problem and, on the other hand, the 
Poisson structure is derived from the classical r-matrix related to the same Rie- 
mann problem, one anticipates some relation between the two matters. However, 
the simplest guess appears to be wrong: dressing transformations do not preserve 
Poisson brackets on the phase space of integrable PDE's. The exact relation is 
more subtle: the loop group itself carries a natural Poisson bracket, again derived 
from the classical r-matrix, and dressing is an example of a Poisson group action 
[18]. The same Poisson structure on the loop group plays an important role in the 
study of difference Lax equations, and also yields a semiclassical approximation in 
the theory of Quantum groups. 

In this paragraph we shall discuss the simplest facts about dressing transforma- 
tions. We start with the motivation of the main definitions. 

Let G be the double loop group consisting of maps 

5 : M X CPi\ ({0} U {oo}) ^ GL (n) 

which are holomorphic with respect to the 2nd argument. The Lie algebra g of G 
consists of maps 

U:Rx CPi\ ({0} U {oo}) ^ qI (n) 
whicli arc lioloniorpliic w ith respect to the 2nd argument. Informally, we may call 

-"-^ There is a diflfcront version of dressing whieh uses the Riemann problem in the half-plane; 
under some additional restrictions, this version is adapted to the study of rapidly decreasing 
solutions. 
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elements of G wave functions. Define the mapping 

G acts on itself by left multiplications; we have 

(8.36) g ■ Ug^ = gU^g'^ + d^g ■ g~^, 

in other words, left multiplication on G induces gauge transformations on the set 
of "potentials" U. Conversely, if ?7 G g, let ipu be the fundamental solution of the 
differential equation on the line 

(8.37) '^ = u{x,z)i,, 

normalized by tpu (0) = Id. The mapping 

■0 : g ^ G : U I— > Vu 

is a right inverse of p. Let ^ C G be the subgroup consisting of maps which do 
not depend on a; e M. 

Lemma 8.30. (i) Q <Z G is the stationary subgroup of the zero potential on the 
line, (ii) Q'^ = ipuG'^u^ CG is the stationary subgroup of U G g. 

Let G+ C G be the subgroup consisting of functions which are regular in 
CPi\{oc} with respect to the second argument, and G_ C G the subgroup of 
functions which are regular in CPi\{0} and satisfy the normalization condition 
g- (oo) = Id. The factorization problem in G consists in representing g G G as 
g = g+gZ^,g± € G±; the first argument of g{x,z) is regarded as a parameter. 

Theorem 8.31. Formula 

(8.38) dr(^^^y)-il) = {tpxy~^ip~'^)_^^ ipx = {ipxy~^t(j~^) _^ tpy 

defines a right group action dr : {Q x Q) x G > G. 

The definition looks rather exotic; in particular, the composition law for dressing 
transformations is not at all obvious. Note that the equality in (8.38) is closely 
related to the fact that -ipxy^^^/j"^ G Q^^' and hence the two factors {^Jxy^^ijj^^) _^ 
define the same gauge transformation of the potential U^p. To check the compo- 
sition law we shall give a geometric interpretation of (8.38). To avoid lengthy 
notation we shall use a model example. 

Let if be a group admitting a factorization into product of its subgroups K± ; set 
D{K) = KxK. Let C KxKhe the diagonal subgroup, = {{x, x);xe K}, 
and Kj. = x . 
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Lemma 8.32. D{K) = Kr • ; in other words, the factorization problem in 

D{K), 

(8.39) {x, y) = (r?+, 7?-) • (€, ^) , n± & K±, ^ & K, 

is uniquely solvable}'^ 
Indeed, we have 

ri±= {xy~^)^, £, = {xy-^)~^x= {xy-^)~^ y. 

CoroUciry 8.33. The quotient space Kr\D{K) of left coset classes modKr is 
modelled on the diagonal subgroup ; the projection it : D (K) is given by 

Tr{x,y) = {xy~'^)^^x= {xif^)_^y. 



The group D{K) acts on itself by right translations. Consider the commutative 
diagram 

D{K) X D{K) D{K) 



TT X id 
Kr\D{K) X D{K) 



TT 



Kr\D{K) 



K X D{K) — K 

Proposition 8.34. The right action K x D{K) K induced by the identifica- 
tion of C D{K) with the coset space Kr\D{K) is given by 

(8.40) dr{x,y) : k i — > {kxy~^k~^)_^kx = [kxy~^k~^) _ky. 

A comparison of (8.38) and (8.40) explains the mystery around the definition. 

In our model setting we assumed for simplicity that the factorization problem is 
globally solvable. In general this is of course not true; however, under reasonable 
conditions it is solvable on an open dense subset of the big group, and hence the 
diagonal subgroup may be identified with a "big cell" in the qiioticnt space. Thus 
the situation is not much different from the treatment of e.g. fractional linear 
transformations on the line. 

Returning back to (8.38), note that the diagonal subgroup C T) [Q) acts by 
dr {g,g) : 4' this action amounts to a simple change of the normalization 

^'^D(K) is called the double of K; we have already used a similar construction for Lie algebras 
in Section 4. One may notice tliat the construction below uses only factorization in D(K) and 
the subgroups K± need not be complementary in i^. In Section 9 we shall once again encounter 
this construction in the study of Poisson Lie groups. 
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of the wave function and does nor affect the potential U = dxfpilJ~^- On the 
other hand, the subgroup Gr = G+ ^ G- preserves the normalization condition 
tp (0) = Id. Hence we may define an action Gr ^ S ^ S on the space of potentials 
with the help of commutative diagram 

X G — 



dr 



(8.41) 



id X tp 

GrXS 



dr 



g 



Let gM.Jv C g be the subspace of Laurent polynomials with the degrees of pole 
at zero (at infinity) not exceeding M (resp., A'') 



Proposition 8.35. Dressing action on g preserves gM,N- 

Sketch of a proof. Compare two equivalent formulae for dressing which follow 
from (8.38); the first one shows that dressing does not increase the degree of pole 
at 0, the second one, that it does not affect infinity. 

This argument explains the key idea of the dressing method: indeed, the most 
striking property of dressing is the fact that it preserves the structure of poles of 
the Lax operator. A slight refinement of the same argument shows that dressing 
preserves symplectic leaves of the r-bracket in gM,N C g ~ g* (here r is the 
standard r-matrix associated with the factorization problem in G). 



Proposition 8.36. Dressing preserves strong regularity. 

Sketch of a proof. Formal Baker functions at zero and at infinity of the dressed 
operator are given by 

= (%+5zv-'):'*oo. 

Clearly the gauge factors (V'.9+.9l^V'^^)_|_^ j (V'5+5-^V'^^)_^ expand into conver- 
gent series in local parameter around zero and infinity, respectively; hence the 
same is true for the dressed wave functions. 

In applications, dressing is usually applied to trivial, or free. Lax matrices. Let 
us assume that the leading coefficient at infinity is a diagonal matrix with distinct 
eigenvalues. By definition, free Lax operator has the form 

Lfree = ^ ~ •^(-^)' 

where D{z) is a constant diagonal matrix (with coefficients which are polynomial 
in z). Our next assertion shows that the factorization theorem survives for reg- 
ular potentials; moreover, the dynamical flows associated with all Lax equations 



Integrable Systems and Factorization Problems 



51 



(derived from a given Lax operator) correspond to the action of an abelian sub- 
group of the "big" dressing group (essentially, the group of diagonal loops which 
are regular at infinity). 

Proposition 8.37. Assume that L is obtained from Lfree by dressing, L = L^f^^^. 
The integral curve of the Hamiltonian equation of motion with the Hamiltonian 
Ha defined (8.27) which starts at L is given by 

L{t)=g±{t)-^oLog±{t), 

where g±{t) are regarded as multiplication operators on the line and g+{t,x), 
g- {t, x) are the solutions of the factorization problem^^ 

g+{t,x)g_{t,x)~'^ = tpfreeix)expta{z~'^) -g-ipfreeixy^ , 

i>free{x) = expxD{z). 



9. Difference Equations and Poisson Lie Groups 



9.1. Motivation: Zero curvature equations on the lattice 



We have already mentioned that integrable systems which are associated with 
difference operators require a special treatment; in this case the underlying Poisson 
structures are nonlinear, and hence the geometric setting we considered so far, 
based on the use of the Lie-Poisson brackets, must be generalized. Nonlinear 
equations associated with a finite difference operator may be regarded as lattice 
analogues of zero curvature equations. They are usually written in the form 

(9.1) ^ = LmMm+l - MmLm, TU € Z. 

Equation (9.1) is the compatibility condition for the linear system 

(9.2) i^m+l = Lmtprn, 

This system of equations is covariant under the gauge transformations of the form 

(9.3) Ipn gmtpm, i"^ gm+lLrnQm ^ ^rn ^ gmMmg^ + ^tgm ■ 9^ ■ 



^^This problem is nontrivial, since in order to get a nonzero Hamiltonian Ha, the constant 
diagonal matrix a (^~^) must have pole at zero. 
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The use of difference operators associated with a 1-dimensional lattice is par- 
ticularly well-suited for the study of "mTilti-particlc" problems. Let us assume 
that "potentials" are periodic, im+JV = L^; the period N may be interpreted 
as the number of copies of an elementary system. In this way we get families 
of Hamiltonians which remain intcgrablc for all A''. The phase spaces for such 
systems are direct products of "one-particle" phase spaces. 

It is natural to suppose that the dynamics associated with difference Lax equa- 
tions develops on submanifolds of a matrix Lie group G (or of a loop group, if there 
is a spectral parameter), rather than on Lie algebras or their duals. As before, 
we are looking for a geometric theory which should simultaneously account for 
the Poisson structure of the phase space, the origin of conservation laws, and the 
reduction of dynamics to factorization problems. An extension of the geometric 
scheme described in Section 3 to lattice systems is based on the theory of Poisson 
Lie groups introduced by Drinfeld [20] following the pioneering work of Sklyanin 
[19].^^ Unlike the Lie-Poisson brackets discussed before, this new class of Poisson 
brackets was virtually unknown in geometry. 

Very briefly, the motivation for the formal definitions we are going to discuss 
is as follows. As in the continuous case, the natural Hamiltonians associated 
with the zero curvature equations should be gauge invariant. Let us assume that 
Lm+N = Lm- Consider the monodromy mapping which assigns to the set of local 
Lax matrices their ordered product. 



T : —>■ G : {Lq, Lat-i) ' — > Tl =JJ L^. 

k 



A version of the Floquet theorem asserts that two difference operators with peri- 
odic coefficients are gauge equivalent if and only if their monodromy matrices are 
conjugate. Thus one expects the integrals of motion of equation (9.1) to be of the 
form 

hk=tvTl 

This will hold if the monodromy itself satisfies a Lax equation, 
(9.4) ^ = l^--*'l 

In more formal terms, let Ft : G^ — > G^ be the dynamical flow associated 
with (9.1) and Ft : G — >G the corresponding flow associated with (9.4); then the 



^^Typical dynamical systems of this kind are the classical analogues of lattice models in quan- 
tum statistics, although some of the systems which we mentioned earlier (e.g., Toda lattices, 
certain tops, etc.) also admit difference Lax representations. 

■''^Tliis latter paper was in turn motivated by the Quantum Inverse Scattering Method devel- 
oped by Faddeev, Takhtajan and Sklyanin (cf. [22]) and by the work of Baxter on Quantum 
Statistical Mechanics (cf. [21]). 
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Tl 



Tl 



(9.5) 



G 



Ft 



G 



We want to equip our phase spaces with Poisson structures which are compatible 
with all mappings in this diagram. Moreover, we would like to keep to our geomet- 
ric picture suggested by theorem 3.1; this means that we must find two Poisson 
brackets in each space, so that 

(i) Spectral invariants of the monodromy are Casimir functions for the first 
structure. 

(ii) They are in involution with respect to the second one and generate difference 
Lax equations (respectively, ordinary Lax equations for the monodromy). 

(iii) The flows Ft, Ft preserve intersections of symplectic leaves for the two brack- 
ets. 

(iv) Vertical arrows in the diagram (9.5) are Poisson mappings with respect to 
both structures. 

(v) Finally, the equations of motion (both upstairs and downstairs) should reduce 
to a factorization problem. 

It is remarkable that all these conditions may be satisfied with the help of a 
single ingredient, the classical r-matrix, the same one which is responsible for the 
factorization problem. As compared to the previous case (that of Lie algebras), we 
need only one extra property (which actually is satisfied in most of the examples 
we considered beforehand): the Lie algebra Q of our Lie group G must he equipped 
with an invariant inner product and the r-matrix r G EndQ must be a skew- 
symmetric operator. The construction which provides Poisson brackets satisfying 
all these conditions is rather nontrivial (in fact, an important message is that this 
is possible at all!); it may be naturally divided into two separate problems: 

1. Given an r-matrix, find the brackets on G = G^ and on G which have 
spectral invariants of the monodromy as their Casimir functions. 

2. Find a Poisson bracket on G = G^ which yields zero curvature equations 
(9.1) as the equations of motion and assures that the monodromy map is 

compatible with the Poisson brackets. 

The key point in both questions is that the r-matrix is fixed in advance and we 
must arrange the brackets with its help (otherwise, there are too many options 
and the problem is not well posed!). 
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The second question is better known than the first one; in fact, it is this ques- 
tion that has led to the theory of Poisson groups. The key step is the following 
simplifying assumption: 

• Dynamical variables associated with different factors in G = G x ... x G 

N 

commute with each other. 

By induction, the monodromy T : G^ — > G is a Poisson mapping if the product 
map m: G X G ^ G preserves Poisson brackets. 

Definition 9.1. Poisson bracket on a Lie group G satisfying the property above 
is called multiplicative; a Lie group equipped with multiplicative bracket is called 
a Poisson Lie group. 

Let us explain this condition in a more explicit way. Let ip,ip & C°°{G); put ^{x,y) = 
fixy), '${x,y) = Pixy), * G C°°{G X G). In order to compute the Poisson bracket {3>, ^} we 
regard them as functions of two variables, that is, we compute derivatives of <E>, with respect 
to X for fixed y and with respect to y for fixed x and take the sum of these two terms; on the 
other hand we may compute the bracket {ip, ip} for functions of one variable z E G and then 
insert z = xy. Multiplicativity means that the two results coincide. 

9.2. Key example: Sklyanin bracket 

Let us assume that the Lie algebra of G carries an invariant inner product and 

r £ End g is skew and satisfies the modified Yang-Baxter identity. For tp € C°° (G) 
let V(p, VV G be its left and right gradients defined by 

{Vv{x),X) = (^^) ^^^{e^^'-x), 

(VV(a;),^) = _^v(a;-e^^),xe0. 

Proposition 9.2. The bracket on G 

(9.6) ^} = ^ ((r (V<^) , VV') - (r (W) , V») 

is multiplicative and satisfies the Jacobi identity. 

^^Formula (9-6) seems to be not the simplest bracket which can be arranged using an antisym- 
metric operator: why not taJje {ifjip}^ (x) = (r ( V<^) , V?/>) , or {(^, (x) = (r (VV) , V^) ? 
The reason is this: when r satisfies the modified Yang-Baxter identity, neither of these brackets 
satisfies Jacobi. However, the obstructions cancel each other when wc take the difference, or the 
sum of the two (and precisely in these two cases)! We shall return to this question in Section 
9.4 below. 
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Remcirk 9.3. Formula (9.6) coincides with (8.21), which we deduced from the 
study of the monodromy map in the continuous case. 

This is of course not a coincidence. To explain why the Poisson bracket for the monodromy on 
the circle should be multiplicative let us consider the auxiliary problems (8.6) with potentials L 
consisting of two separate patches, so that supp L is the union of two disjoint intervals, supp L = 
I'Ul" . Let us denote by 25// , 25/// the set of all potentials supported on /', /", respectively; then 
25//, 25/// C 25 are Poisson submanifolds, and moreover, 25/'^/// = 0// X 0/", again as Poisson 
manifolds (which means that L' G &ji,L' G 25// may be treated as independent variables with 
respect to our Poisson structure). Clearly, for L = L' + L" , L' G (3ji,L" G 25/// we have 
M/, = Mi^iiM]^! and the Poisson bracket for the monodromy matrix may be computed in two 
different ways: either by computing the Poisson brackets for the monodromy matrices M/^//, Mj^i 
regarded as independent variables, or alternatively by computing directly the monodromy M/^ 
for the potential L = L' + L" decomposed into two separate patches. The two results should of 
course coincide, and this means precisely that the Poisson bracket for the monodromies should 
be multiplicative. 

Note that Vtp {x) = Adx ■ VV (a^)> or, in the matrix case, Vip {x) = x ■ VV (a;) • 
x~^, so we may rewrite (9.6) using only left gradients: 

{^, V} {x) = tr (r?, {x) ■ {Vip A W)) , 

where we set 

rjr [x) = — {r — Adx~^ or o Adx^ 

and identify Vip A VV' G A g with an antisymmetric linear operator on fl, using 
the inner product. The function rjr : Q ^ Endfl satisfies the following remarkable 
functional equation:^° 

(9.7) rjr (xy) = rjr (x) + Adx~^ o rjr (y) o Adx, 

One may check that this functional equation is exactly equivalent to the multi- 
plicativity of (9.6). 

Assume that G = GL (n) is a matrix group. Let us consider "tautological'" 
functions (pij on G which assign to a matrix i e G its matrix coefficients, (L) = 
Lij-, clearly, the ring of polynomials C [(f>ij] is dense in C°° (G), and the Poisson 
bracket on G is complcitcly specificid by its values on the "generators" (pij. Let us 
identify r € EndBt(n) with an element of 5t(n) (8) gl{n) ~ Mat (n^) . 

Proposition 9.4. The Poisson bracket (9.6) of matrix coefficients is given by 

(9.8) (L) = [r,L® L\^.^ . 

Pit this point we use the crucial property of (8.11): the Poisson operator is a multiplication 
operator. 

^°This condition is usually expressed by saying that rir is a 1-cocycle (in fact, a coboundary) 
on G. In these lectures we shall not use this language. 
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The commutator in the r.h.s is computed in Mat (^^^)■ Usually people do not 
distinguish <pij and its values and write this formula (with suppressed matrix 
indices!) as 

(9.9) {L f L} = [r,L^L]. 

Formula (9.9) has served as the original definition of the Sklyanin bracket. Note 
that the r.h.s. in (9.9) is a quadratic expression in matrix coefficients (this is to 
be compared with the Lie-Poisson bracket which is linear). 
Let us note some important properties of the Sklyanin bracket. 

(i) The bracket is identically zero at the unit element of the group. (Indeed at 
X = e right and left gradients coincide). 

(ii) Linearizing the bracket at the origin of group, we get the structure of a 
Lie algebra in the cotangent space T*G — Q* : li Q* ,X ^ q, choose 
fjipG C°° (G) in such a way that Viy? (e) = 5, (e) = Q and set 

(the second formula checks that the bracket C]* does not really depend on 
the choice of ■0 and so is well defined). 

Proposition 9.5. The bracket C]* coincides with the r-bracket (up to the iden- 
tification of and g* induced by the inner product?^ : 

(9.10) [^,Cl = liH,C] + [trC]). 

Up to dualization, (9.10) coincides with (3.1) which was our starting point in 
Section 3. In the present setting we get some extra properties which follow from 
multiplicativity of the bracket. Set 

Explicitly, we get 

Sr {X) = adX or — r o adX. 
Proposition 9.6. (i) We have 

(9.11) tr(5,(X)o(eAC)) = ([e,C],,^). 

Formula (9.10) explains the choice of normalization in (9.6): we wanted to get the same 
thing as in (3.1). 
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(ii) The mapping 6r : Q ^ Endg satisfies the functional equation 
(9.12) 6r {[X, Y]) = [adX, 6r {¥)] - [adY, 6r {X)] . 

Equation (9.11) shows that 6r is the dual of the commutator map g* Afl* —>■ g*.^^ 

Definition 9.7. A pair (0,0*) is called a Lie bialgebra if (i) g and g* are set in 
duality as linear spaces, (ii) both g and g* are Lie algebras, (iii) the dual of the 
commutator map [, ]^ : g* A g* — > g* satisfies the functional equation (9.11). 

Remark 9.8. One can show that (iii) implies that in, the dual way, the mapping 
: g* — > g* A g* which dualizes the commutator [, ] : g A g — > g is a 1-cocycle on 
g*, and so this definition is symmetric. 

It is instructive to compare the definitions of Lie bialgebras and of the double Lie 
algebras introduced in Section 3. These definitions are different and use different 
notions of the classical r-matrix. In the case of double Lie algebras there are 
two Lie brackets on the same underlying linear space; the classical r-matrix is a 
linear operator r G Endg; in the case of Lie bialgebras there are two Lie brackets 
which are defined on dual spaces g and g*. The motivation for these definitions 
are very much different as well: as we saw, double Lie algebras provide a natural 
setting for the Involutivity theorem (theorem 3.1); Lie bialgebras naturally arise 
in the study of multiplicative Poisson brackets on Lie groups. Proposition 9.5 
specifics the setting in which these two notions merge together: we must assume 
that g carries an invariant inner product which allows to identify g and g* and 
that r € Endg is skew. One more natural condition is the modified Yang-Baxter 
equation (which assures that there is an underlying factorization problem) . When 
all three conditions are satisfied, we say that (g, g*) is a factorizable Lie bialgebra. 
Factorizable Lie bialgebras and the associated Poisson Lie groups provide a natural 
environment for all applications to lattice integrable systems. 

Before turning to lattice zero curvature equations let us discuss ordinary Lax 
equations on Lie groups. Here is a version of the factorization theorem (theorem 
3.3) which applies in this setting. Let G be a matrix Lie group; we assume that 
the Poisson bracket on G is given by (9.6) and that its tangent Lie bialgebra 
(g,g*) is factorizable. Let /(G) G C°° (G) be the algebra of central functions on 
G {(f £ C°° (G) is central if ip {xy) = (fi (yx) for all x,y £ G). 

Theorem 9.9. (i) All central functions are in involution with respect to the 
Sklyanin bracket (9.6). (ii) Hamiltonian equation on G with Hamiltonian h G I (G) 



'In formal terms, equation (9.12) means that Sr is a 1-cocycle on g (with values in Endg). 
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may be written in Lax forrn^^ 

(9.13) ^ = LM± - M±L, 

at 

where M± = r± (V/i (L)). (Hi) The integral curve L {t) of (9.13) with L (0) = Lq 

is given by 

(9.14) Lit)=g±it)-^Log±it), 

where g+ (t) , g- {t) are the solutions of the factorization problem in G 

(9.15) g+{t)g-{t)-^ =exptVh{Lo) 
associated with r. 

As before, the direct proof of theorem 9.9 is easy: to check that (9.14) is an 
integral curve of (9.13) just compute the derivative of the r.h.s in (9.14). As in 
Section 3.3, there exists also a geometric proof which explains the background 
machinery. Below, we shall briefly outline the corresponding construction. 

9.3. Duality for Poisson Lie groups 

As already noted. Lie bialgebras possess a remarkable symmetry: if (0,0*) is 
a Lie bialgebra, the same is true for (0*,0). Hence the dual group G* (which 
corresponds to G*) also carries a multiplicative Poisson bracket. In the case of 
factorizable Poisson groups this dual bracket may be pushed forward to G by 
means of the factorization map. Thus we get two brackets on G which fit into our 
geometric treatment of Lax equations. The best way to understand this duality 
is to notice that both G and G* are Poisson subgroups of a bigger Poisson group, 
the double ofG. 

Let (0, 0*)be a Lie bialgebra; the linear space = ® 0* carries a natural inner 
product 

(9.16) {{{X, F) , {X', F'))) = {F, X') + {F', X) . 
The following key theorem was discovered by Drinfeld. 

Theorem 9.10. There exists a unique structure of the Lie algebra on D such that: 
(i) 0,0* CO are Lie subalgebras. (ii) The inner product (9.16) is invariant. 

Corollciry 9.11. Let Pq^P^- he the projection operators onto 0,0* in the decom- 
position = © 0*. Set rr, — Pg — Pg*; then defines on 5 the structure of a 
factorizable Lie bialgebra. 

^^In this formula the velocity vector dL/dt belongs to the tangent space Ti,G; in order to be 
more accurate, we may rewrite (9.13) as an equality in the Lie algebra: 

L-^dL/dt = M±- AdL-^ ■ M± 



Integrable Systems and Factorization Problems 



59 



The pair (0,0*) is called the Drinfeld double of When the initial Lie 

bialgcbra (0,0*) is itself factorizablc, the description of the double is very simple. 
Consider the Lie algebra = g © g (direct sum of two copies of g) and equip it 
with the inner product 



where (, ) is the invariant inner product on q. 

Proposition 9.12. The double of a factorizable Lie algebra is canonically iso- 
morphic to d = © 0. 

Sketch of a proof. We have already seen in Section 4 that there arc two natural 
homomorphisms r± : 0* ^ given by (4.5); their combination yields an embedding 
0* C © 0. Let 0"^ C © be the diagonal subalgebra, 0*^ = {{X,X)-, X G 0}. 
As discussed in Section 4, = 0'^+0*; it is easy to check that the skew symmetry 
of r and the choice of the inner product in imply that 0* and 0"^ are isotropic 
with respect to the inner product (9.17): this is equivalent to the skew symmetry 
of rr, = Pg — Pg. . In matrix notation, rj, G End (0 © 0) is given by a 2 x 2 block 
matrix: 



The Lie group which corresponds to is D (G) = G xG. Let G^, G* be its sub- 
groups which correspond to 0'',0*. Clearly, G^ C D (G) is the diagonal subgroup. 
As in (8.39), we may associate with the r-matrix rj, a factorization problem in 
D (G). Let us assume for simplicity that it is globally solvable, i.e. D (G) c^iG-G* . 

Proposition 9.13. Let us equip D [G) with the Sklyanin bracket associated with 
r-t,- Then G^,G* C D (G) are Poisson subgroups (i.e. they are Poisson submani- 
folds and the induced Poisson structure is multiplicative) .'^'^ 

The bracket induced on G^ coincides with the original Sklyanin bracket associ- 
ated with r; the Poisson bracket on the dual group G* is described in the following 
way. Consider the mapping m : D{G) G : {x,y) ^ xy~^; its restriction to 
G* C D (G) is a diffeomorphism. 

Proposition 9.14. The Poisson bracket on G induced by m : G* ^ G is given by 



where Vip,Vtp and VVi VV are left and right gradients of ip,tlj. 

Formula (9.19) looks rather complicated; however, the bracket {, }^ is very re- 
markable. 



(9.17) 



{{{X,Y),{X\r))) = {X,X')-{Y,Y') 



(9.18) 




(9.19) 



= l/2(rV^,VV') + l/2(rVV,V» 
- (r+ V^, V» - (r_ W, VV') , 



More precisely, the bracket induced on G* C D (G) has opposite sign, due to the minus sign 
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Proposition 9.15. (i) Sym,plectic leaves o/{,}^ coincide with conjugacy classes 
in G. (ii) Casimir functions of {. }^ are precisely the central functions on G. (Hi) 
The bracket (9.19) vanishes at the unit element e e G; the induced Lie bracket on 
the tangent space coincides with the original Lie bracket on q. 

Thus the bracket (9.19) provides the missmg ingredient of our geometric picture: 
we have got now two Poisson structures on the same underlying manifold G and 
Lax equations preserve intersections of two systems of symplectic leaves. 

The symplectic leaves of the Sklyanin bracket also admit a description in terms 
of the factorization problem. Let us identify G with the quotient space D (G) / G* 
using the factorization D (G) = G^ -G* Let us denote by tt the canonical projection 
D (G) — > D (G) /G*; define the action G* xG ^ G using the commutative diagram 



G* X D{G) 



id X 



G* X D{G)/G* 



m 



D{G) 



D{G)/G* 



G* X G '^^^'"'^ ' G 

(Here m is the group multiplication in D (G) restricted to the subgroup G* C 
D (G) .) By analogy with the definition of dressing transformations, this action is 
called dressing action. 



Proposition 9.16. Symplectic leaves of the Sklyanin bracket in G coincide with 

the orbits of G* in G with respect to the dressing action. 

More explicitly, the dressing prescription is as follows: given x & G, {h^,h—) € G* solve the 
factorization problem in D (G); 

{h+x,h-x) = [x'g+,x'g-^ ,x' e G,{g+,g-) e G*: 

then Dress {h^,h—) ■ x = x' . This immediately yields the following formula in terms of the 
factorization problem in G : 

(9.20) Dress (h-f-, h-) ■ X = h-f-x (^x^'' h^^'h-x^ ^ = h-x (x~^h^^h—X^ , 

where x~^h^^h-X = (^x~^h^^h-X^ ^ ■ (^x~^h^^h-x^ is the factorization in G associated 
with the original r-matrix. 



Exercise 9.17. Check this formula using the definition of dressing action. 

Dressing action may be regarded as a nonlinear analog of the coadjoint repre- 
sentation, as it is clear from the following simple assertion. 



Proposition 9.18. Dressing action leaves the unit e € G invariant; the lin- 
earization of the dressing action in the tangent space T^G c:^ q coincides with the 
coadjoint representation of G* in q** ~ q. 
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In applications it is natural to assume that G is an algebraic loop group consisting 
of matrices whose coefficients are rational functions of z. Orbits of the dressing 
action of G* in this loop group are finite-dimensional, and we get a description of 
phase spaces for Lax equations which is largely parallel to the case of coadjoint 
orbits of Gr discussed in Section 6. 

One more application of the dual Poisson structure described by (9.19) is the 
accurate description of the Poisson properties of the dressing transformations from 
Section 8.6. We keep to the notation introduced in lemma 8.30. Let Q be the loop 
group; its Lie algebra Lq = q[z, z~^] is equipped with the inner product 

(X,y)=Res,=otr(X {z)Y [z)) 

and with the standard r-matrix r = P+ — P_ associated with the Riemann problem. 
The loop group Q equipped with the corresponding Sklyanin bracket becomes a 
Lie-Poisson group. Let Gr — G* — G+ x G- be the dual group equipped with its 
natural Poisson structure (9.19). 

Theorem 9.19. The dressing action described by the commutative diagram (8.41) 
is a Poisson group action. 
(We shall not reproduce the proof here; see [18], [23].) 

9.4. Symplectic double and the free dynamics 

One more important ingredient of the geometric picture outlined in Section 3.3 
is the "big phase space" with "free" dynamical flow. Its counterpart in the present 
setting is provided by the so called symplectic double of G. Let again D (G) = GxG 
be the double of G; the Sklyanin bracket on D (G) is given by 

(9.21) V4 = ((ro V^, V^A)) - ((r, W, V^)) - 

As noticed in footnote 18, the terms with left and right gradients separately do 
not satisfy the Jacobi identity; the obstructions cancel when the two are combined 
together. Explicitly, 

{{^1, ip2}^, ip3}<^ + C.p. = (([V^l, V^2] , V^s)) , 

(9.22) 

{{^1, (^2}^ ^3}'' + C.p. = - (([Wl, VV2] , VVs)) . 

The two terms cancel, since V<^ (x) = x ■ VV {x) ■ x~^ and the inner product is 
^d-invariant. It is important to notice that the crucial minus sign in (9.22) is due 
not to the minus in (9.21), but rather to the fact that the action of a group by right 
translations is its a«ij-representation. (There are no terms of "mixed chirality" in 
the obstruction, since left and right translations commute with each other.) Thus 
we get the following assertion: 
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Proposition 9.20. Let r,r' e EndO be two arbitrary classical r-matrices satisfy- 
ing the modified Yang-Baxter equation; the bracket 

(9.23) = {{rVip, Vw)) + {{r'Vip, VV)) 

satisfies the Jacobi identity. 

Specifically, let us take r = r' = rj,; the resulting Poisson structure is non- 
degenerate (at least if wc assume - as we always do in this Section that the 
factorization problem in D is globally solvable), and hence defines a symplectic 
structure on D (G) . 

Definition 9.21. The manifold D (G) with the Poisson bracket 

(9.24) {^, = {{r^Wip, W)) + {{r^Vip, W)) 

is called the symplectic double of G. We shall denote it D {G)_^ in order to distin- 
guish it from the Drinfeld double equipped with the Sklyanin bracket. 

The bracket (9.24) is not multiplicative, so from the point of view of the theory 
of Poisson groups D (G)_^ is not a group! Instead, we have the following property 
which shows that D (G)_|_ is a principal homogeneous space for D (G) : 

Proposition 9.22. Left and right multiplication in D (G) induce Poisson map- 
pings D{G)xD (G)+ — > D (G)+ , D (G)+ x D (G) — >D (G)_^. 

Proposition 9.22 paves the way to use the reduction technique in our present 
setting. Let us recall the point of view on reduction adopted in Section 3.3: if M 
is symplectic and K x M — > M is a group action, the reduction is the natural 
projection map n : M — > M/K onto the space of i^-orbits in M. The key 
property which we need to get a Poisson bracket on M/ K is this: Poisson bracket 
of two G-invariant functions on M is again G-invariant. Let us discuss briefly 
how can one control this property. For X € 4 let us denote by X G VectM the 
vector field on M generated by the 1-parameter transformation group cxptX. We 
have: 

1^ g G°° (M) is G-invariant Xip^O for all X €t 
When vector fields X e Vect M are Hamiltonian we have simply 

xw,i,} = {x^,^} + y,x^} = o. 

In the case of Poisson group actions vector fields X are no longer Hamiltonian; 

however, the rate of nonconservation of Poisson brackets by these vector fields 
may be characterized very sharply. For G C°° (M) , x G M, let us denote by 
^tp (x) G 6* the linear functional defined by 
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Proposition 9.23. Let us assume that K is a Poisson Lie group with Lie bialgebra 
(6, t*). The mapping K x M — > M is a Poisson mapping if and only if 

X{^, V} - {X^, i;} - Xi;} = {[^^,^^1 , X) . 

When Xip = = for all X e t,^^ = = 0, and hence X {ip,ip} = 0, 
which assures the possibility of reduction. In a more general way, let us say that 
a subgroup H d K ]s admissible if Xip = Xip = for all X G ^ X {(p, ip} = 0. 

Proposition 9.24. H C K is admissible <^ [)-'- C t* is a Lie subalgebra. 

Sketch of a proof. Xcp = Xip = for all X G f) imphes that ^<^,^^ G f)^. 
When f)-"- c 6* is a Lie subalgebra, {[^,p,^^]^ ,X) = for all X G f) and hence 

xW,iP} = o. 

As a first example of reduction, let us derive from (9.24) the dual Poisson bracket 
(9.19) on G. 

Proposition 9.25. Consider the action G x D{G)^ — > D{G)^ : h : {x, y) i — > 
{hx,hy). This action is admissible, and the projection map p : D{G)_^ — > 
G : {x,y) I — > x~^y is constant on its orbits. The quotient Poisson bracket on 
D {G)_^ /Gc^G coincides with (9.19). 

Here is a simple check of this assertion: Let ip S C°° (G); set ^{x,y) = (^x^^yj . Clearly, 
^ is left C'-invariant; hence its left gradient V"!" is in (fl*) C 0; but the inner product in 5 
is so chosen that fl'' C is a maximal isotropic subspace, i.e., (fl*) = 2^ ■ Let us compute the 
Poisson bracket $2}_|_ for two functions of such type. Since V$i, V$2 € 0* and rjigs = id, 
we have 

((raV$i, V$2» = ((V$i, V*2)> = 0. 

On the other hand^^ , 

substituting this expression into ((rj V'$i, V'$2)) and using (9.18), we get (9.19). 

We can now state the nonlinear version of reduction theorem described in Section 
3.3. Let (p £ I (G) be a central function. Define the Hamiltonian h^ on D (G)_|_ 
by hy, {x,y) = ip (x-'^y). 

Lemma 9.26. (On free dynamics) The Hamiltonian flow on D{G)_^_ defined by 
h^ is given by 

(9.25) Ft : {x,y) ^ (xe*^, ?ye*^) , X = Vip{x-^y) . 



^^Gradients in the l.h.s are computed with respect to two variables x,y; in the r.h.s they are 
computed with respect to a single variable. 
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Theorem 9.27. (i) The Hamiltonian is invariant with respect to the group 
G* which is acting on D via 

{h+,h-) : {x,y) i — > {h+xhZ^ , h+yhZ^) 

(a) G* is an admissible subgroup in D (G) x D (G) (which acts on D (G)_,_ by left 
and right translations). (Hi) The mapping 

TT : D (G)+ -^G:{x,y)^ y-'xy. 

is constant on G* -orbits in D{G)j^ and allows to identify the quotient space 
D{G)^/G* with the subgroup G = {{x,e);xG G} C D{G)_^_.^^ (iv) The quo- 
tient Poisson structure on D{G)_^/G* ~ G coincides with the Sklyanin bracket, 
(v) The quotient flow Ft on G is given by Ft : x i — > g± it)~^ xg± (t), where 
g+{t),g-{t) solve the factorization problem exptV(fi {x) = g+ {t) g- {t)~ , and 
satisfies the Lax equation (9.13). 

9.5. Lattice zero curvature equations and the twisted double 



So far our geometric construction is restricted to ordinary Lax equations on 
a single copy of G. In order to put lattice zero curvature equations into our 
framework wc need one more effort. First of all, let us state the factorization 
theorem which applies in this case (cf. Section 9.1). Let again 

T : G^ — > G : (Lq, L^-i) ' — > 7l Lk 

k 

be the monodromy map; choose (fi € I (G) and set Hip — (p o T. Wc define the 
"wave function" tprn associated with the auxiliary linear problem (9.2) by 

0<fe<m-l 

The Poisson structure on G^ is defined as the direct product of Sklyanin brackets 

on each factor. As usual, we assume that our basic r-matrix is skew and satisfies 
the modified Yang-Baxter equation, so that (fl,fl*) is a factorizable Lie bialgebra. 

Theorem 9.28. (i) The Hamiltonian equation of motion on G^ with Hamiltonian 
Hip may be written as 

(9.26) ^-^=L^M^^,-M^L^, 
^® As ijsual, for p g G we denote by g+ , p_ the solutions of the factorization problem g = g+gZ^ . 
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where 




where g+ (t) , g, 



{t) are the solutions of the factorization problem 



(9.29) 



(t)-^=^^-exptV9J(rio).«) 



-1 



As usual, a direct proof of theorem 9.28 is easy. One point worth to be mentioned 
is the computation of the gradients of the Hamiltonian: for that end we may use the 
"variation of constants" in the auxiliary linear problem (9.2), similarly to the case 
of differential operators on the circle discussed in Section 8.3. A geometric deriva- 
tion is not so straightforward. Let us introduce the following notation in order to 
simplify the bulky formulae. Wc set G = G^,(S = ©^0,L = (Lo, Ln-i) e 
G. Let r be the automorphism of G induced by cyclic permutation of indices, 
{Lq, ...,Li^-\Y = {Li,L2, ...,Lj\f-i,Lo); we denote the corresponding automor- 
phism of (S by the same letter.^'' Equations (9.26), (9.28), (9.29) may be rewritten 
as 



(in the last expression exptVip is embedded into G = diagonally). The 
important feature of these formulae is the presence of the twisting automorphism 
T. It plays the key role in the definition of two more objects: (1) The second 
Poisson bracket on G which has got the spectral invariants of the monodromy as its 
Casimirs and (2) The twisted Poisson structure on the "big double" D (G) = G x G 
which is responsible for the free dynamics. We shall start with the latter question. 

Let r S End g be the classical r-matrix associated with the factorization problem 
in a single copy of G. The corresponding r-matrix in the "big" algebra g is the 
direct sum of r's acting in each copy of q, 



Note that roT = Tor. Ina similar way, the r-matrix of the "big double" S = ©©25 
is 



(9.30) 



f = LMJ - M±L, 
L(t) = g± (t)-^L°g± (t)^ 



where 



g+ (t) g^ (t)-^ = ^° cxp tV^. (Tlo) (^°) 



-1 



^''The twisting automorphism r plays the role which is similar to that of the derivation dx for 
ordinary zero curvature equations; we can say that its use allows to reproduce for lattice systems 
the effects of central extension of loop algebras discussed in Section 8.1. 
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in block notation we have 




(cf. (9.18)). Let us define the automorphism r e Aut (G x G) by (x, y)'^ = 
(x, y"^); the corresponding automorphism of 2) = 6 © 25 is again denoted by the 
same letter. Put rj, = r o o r"-^; in block notation we have 

r 2r+ o T^-*- \ 

2r o r_ — r J ' 

Let us define the twisted Poisson structure on D (G) by 

(9.31) V'}, = ((rs V<^, VV')) + ((rj W, V») . 

The Jacobi identity for (9.31) follows from proposition 9.20. Our next assertion is 
the twisted version of proposition 9.25. 

Proposition 9.29. The action G x D (G)^ — > D (G)^ : h : (x,y) i — > (hx.hy) 
is admissible; the projection map p : D (G)^ — > G : (x,y) i — > x~-^y is constant 
on its orbits. The quotient Poisson bracket on D (G)_|_ /G ~ G is given by 

M (v.«-((.-i(!^v)-(-v'«)))^ 

The proof is parallel to that of proposition 9.25. The properties of the quotient 
bracket are also quite remarkable. (Note that unlike the Sklyanin bracket on 
G = G X ... X G the bracket (9.32) is non-local, due to presence of the twist r). 
Let us consider the monodromy map T : G — > G. Assume that G carries the 
Poisson bracket (9.32) and G is equipped with the bracket (9.19). 

Proposition 9.30. (i) The monodromy T : G — > G is a Poisson mapping, (ii) 
The spectral invariants of the monodromy are the Casimirs of the quotient bracket. 
(Hi) Define the gauge action G x G — > G : x : g i — » xg (x"^ )""*■; when N is odd, 
its orbits coincide with sym,plectic leaves of the quotient bracket. 

In brief, we see that the bracket (9.32) provides us with the missing ingredient 
for our geometric picture: it is the sought for second bracket which we need for a 
geometric treatment of lattice zero curvature equations. 

Lemma 9.31. (on free dynamics) Let (p G I{G); set h^ = (p o T e (G), 
= h^op & C°° (D (G)) . The integral curves of the Hamiltonian Hy, in D (G)^ 
are given by 




(9.33) (x (t) , y (t)) = (a;oe*^, yoe*^) ,X = Vv (T(x- V)) 
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Theorem 9.32. Set G* = (G*) and consider the action G* x D (G)^ — > 
D (G)^ given by 

(h+,h_) : (x,y) ^ (h+xhiSh+y (h^)-') . 

This action is admissible; the quotient space D (G)^ /G* may be identified with 
G by means of the map 

^ : D (G) ^ G : (x, y) ^ y;'xyL"' 

which is constant on the orbits ofG*, and the quotient Poisson structure coincides 
with the Sklyanin bracket. The flow (9.33) admits reduction with respect to this 
action; the quotient flow is coincides with (9.30). 

Theorem 9.32 fills the last gap in our geometric picture and shows that the 
qualitative behaviour of difference zero curvature equations is the same as in the 
case of linear phase spaces. 
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